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Abstract 

In this paper we give a review of the most general approach to de- 
scription of reference frames, the monad formalism. This approach 
is explicitly general covariant at each its step, permitting use of ab- 
stract representation of tensor quantities; it is applicable also to spe- 
cial relativity theory when non-inertial effects are considered in its 
context; moreover, it involves no hypotheses whatsoever thus being 
a completely natural one. For the sake of the reader's convenience, 
a synopsis of tensor calculus in pseudo-Riemannian space-time pre- 
cedes discussion of the subject, containing many expressions rarely 
encountered in the literature but essentially facilitating the consid- 
eration. We give also a brief comparison of the monad formalism 
with the other approaches to description of reference frames in gen- 
eral relativity. In three special chapters we consider applications of 
the monad formalism to general relativistic mechanics, electromag- 
netic and gravitational fields theory. Alongside of the general theory, 
which includes the monad representation of basic equations of mo- 
tion of (charged) particles and of fields, several concrete solutions are 
provided for better understanding of the physical role and practical 
application of reference frames. These examples lead often to exotic 
consequences, but this only clarifies the situation (e.g., cases when 
a rotating electrically charged fluid does not produce any electric 
field in the co-moving reference frame, or when magnetic charges of 
a kinematic nature arise in a rotating frame) . The topic is discussed 
of cases when it is unnecessary to introduce a reference frame, and 
when such an introduction is essential, depending on the problems 
under consideration. Special attention is dedicated to analogy be- 
tween gravitation and electromagnetism, closely related to the drag- 
ging phenomenon. An analogy is also disclosed between mechanics 
and electromagnetic field theory in a sense that in the Maxwell equa- 
tions in rotating reference frames appear terms of the same nature 
as that of the Coriolis and centrifugal forces. 
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Preface 



During many years I have felt a necessity to have at hand a concise, but 
complete enough, account of the theory of reference frames and its ap- 
plications to the field theory. There are, of course, such publications as 
the books by Yu.S. Vladimirov, by V.I. Rodichev, by O.S. Ivanitskaya, 
and by A.L. Zel'manov and V.G. Agakov (all of them only in Russian), 
containing expositions of such ideas, and miscellaneous journal papers. 
They were however written using languages which differed drastically 
from the traditional guidelines closely related to intuition-nourishing 
symbolics, such as time-derivative, spatial divergence and curl. It was 
obviously indispensable to translate ideas of the reference frames the- 
ory into this more physical language not losing at the same time the 
rigour of our theory — thus, without restricting to any kind of approx- 
imations (which is done quite often in such decent branches of science 
as the relativistic celestial mechanics where concrete applications are 
otherwise, of course, unattainable). 
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Chapter 1 
Introduction 



1.1 A general characterization of the sub- 
ject 

The concept of reference frame was introduced in physics at the early 
stage when its formalization had just begun and before introduction of 
system of coordinates and equations of motion. When Galileo Galilei 
has first spoken about the principle of relativity of motion (see Schmut- 
zer and Schiitz (1983)), he already used in fact this concept. The more 
used it Newton (1962) having spoken of the law of inertia (the First 
Law) and of the concepts of time and space. The concept of non-inertial 
reference frame was also introduced and developed in classical mechan- 
ics, together with the corresponding effects, the centrifugal and Coriolis 
forces (see our discussion of the subject in (Vladimirov, Mitskievich and 
Horsky 1987)). The idea of reference frames has played a fundamental 
role in physics; Galilei connected them closely with the thought experi- 
ment method (Gedankenexperiment, so important in quantum theory), 
and Newton who considered a rotating bucket full of water has laid the 
foundation-stone for Mach's principle. In this paper, we shall consider 
mainly the problem of description and further applications of refer- 
ence frames in relativistic physics. The context of general relativity is 
essential since it has inevitably to be used also in the Minkowski space- 
time if non-Cartesian coordinates are used. This approach is crucial 
for interpretation and evaluation of physical effects predicted by four- 
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dimensional theories formulated irrespective to any concrete choice of 
a reference frame, as well as for determination of physical observables 
already in classical (non-quantum) theory i.e. also for interpretational 
purposes. 

We shall speak primarily of general relativistic approach to reference 
frames, for which a concise review of ideas, methods, and applicational 
results will be given. As general relativity we understand all relativis- 
tic physics taking account of gravitation, i.e. the space-time curvature. 
Reference frames are considered in the scope of classical (non-quantum) 
physics, i.e. in the assumption that the observation and measurement 
procedures do not influence the behaviour of the objects being observed, 
in particular, they do not disturb the space-time geometry. This means 
that the reference body is a test one, while one and the same physical 
system of real bodies and fields may be described simultaneously and 
equally well in any of the infinite set of alternative reference frames. 
Since the reference body points are test particles, they may be acceler- 
ated arbitrarily without application of any finite forces (hence, of any 
fields), — their masses are infinitesimal too. The only limitation on 
their motion is due to the relativistic causality principle: since a refer- 
ence body represents an idealization of sets of measuring devices and 
local observers, its points' world lines should be time-like (lie inside the 
light cone). It is not always possible to describe motion of a system of 
bodies in a co-moving reference frame without involving singularities 
(at least, the bodies should not mutually collide). Therefore we shall 
assume that motion of a reference body is described by a congruence of 
time-like world lines. This leads to a very fruitful employment of ideas 
of relativistic hydrodynamics in the theory of reference frames; thus 
a reference body may accelerate, rotate and undergo variable defor- 
mation, which completely characterize the reference frame properties. 
This mathematical description should be consistent with the assumed 
space-time geometry. The resulting description of arbitrary reference 
frames (including non-inertial ones) and the corresponding physical ef- 
fects, has explicitly general covariant form. 

Dealing with the reference frame concept comes sometimes across 
certain psychological obstacles when one either confuses reference fra- 
me with system of coordinates (non-covariant approaches), or connects 
it to any four- dimensional basis whatsoever. The both extremities, in 
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the author's opinion, come together in the sense of gross vulgarization 
of this merely simple notion; what is much worse, these approaches 
are widely practised, which involves frequently a disregard and even 
complete negligence of reference frames in the most parts of physics. 
We shall revisit this problem in section 2.1 that precedes the description 
of the reference frames calculus. 

It is remarkable that all algebraic relations in electrodynamics are 
the same in the standard three-dimensional form and in the formalism 
of arbitrary reference frames in any gravitational fields. But what is 
more important, the dynamical (differential) laws of electromagnetic 
field theory in non-inertial frames differ from their standard three- 
dimensional Maxwellian form by additional terms which are naturally 
interpreted analogously to the forces of inertia of the classical mechan- 
ics. Moreover, even in the flat Minkowski space-time, non-inertial ref- 
erence frames involve inevitably and rigorously a curved physical three- 
space which in the presence of a rotation is also non-holonomic in the 
sense that there exist only local three-space elements which do not form 
global three-hypersurfaces. Thus the expressions of the forces of inertia 
and their field theoretical generalizations either do not change at all, or 
change inessentially when one turns from general to special relativity 
(only the four-curvature goes). For example, Maxwell's equations re- 
tain in a non-inertial frame in special relativity exactly the same monad 
form they have in general ralativity. Under a monad, we understand a 
time-like unit vector field whose integral lines depict worldlines of the 
"particles" forming a test body of reference (see below). But we have 
mentioned the role of reference frames in interpreting theoretical results 
and evaluating effects the theory predicts, not in a casual remark: the 
very solutions of the particles and fields dynamics equations are in gen- 
eral much easier obtained from their standard four-dimensional form 
using simplifications due to a choice of coordinates which corresponds 
to the symmetries inherent in the problem under consideration. The 
use of reference frames comes only later. 

Some remarks about the structure of this paper. In section 1.2 we 
consider mathematical (mainly geometrical) definitions and relations. 
This is necessary since the notations differ substantially from one pub- 
lication to another, and moreover it is necessary to have at hand some 
exotic formulae. This paper may also be of use for colleagues working 
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in other areas than gravitation theory, and a bit of self-containedness 
would help them to get into the subject. 

In chapter 2, after the already mentioned section 2.1, we give in 2.2 
the basic relations of the monad formalism algebra, including purely 
spatial but nevertheless generally covariant (in the four-dimensional 
sense) operations of the scalar and vector products. In further sections 
of the chapter 2, the differential relations of the monad formalism are 
given: on the one hand, these are important characteristics common 
with those of mechanics of continuous media; on the other hand, differ- 
ential operations as well as the corresponding identities are considered. 

Chapters 3 and 4 deal not only with gravitation, but essentially with 
electrodynamics too. The first of them is dealing with the relativistic 
mechanics of electrically charged particles, the second with dynamics of 
the very electromagnetic field. It is worth specially mentioning section 
3.3 in which the equations of motion of a test massive charge are consid- 
ered from the viewpoint of an arbitrary reference frame, 3.4 where an 
example is given of calculation of the frequency shift effect in two char- 
acteristic cases, and 3.5 giving an opposite example, that of an objec- 
tively determined Killing reference frame which enables a standard def- 
inition of the dragging phenomenon. Maxwell's equations are written 
down in an arbitrary reference frame in section 4.3 where non-inertial 
terms are also considered. In section 4.4 the reference frame formalism 
is applied to physical interpretation of an exact solution of the system 
of Einstein-Maxwell equations where a non-test charged fluid does not 
produce in its comoving frame any electric field whatsoever. In 4.5 
kinematic non-inertial magnetic charge density is shown to appear in 
the space-time of an Einstein-Maxwell pp-wave in a very natural (but 
rotating) frame of reference. In section 5.2 we consider 3+1-splitting 
of Einstein's equations in the monad formalism and trace some analo- 
gies and dissimilarities between the electromagnetic and gravitational 
interactions. Then, in sections 5.3 and 5.4 the next stage of analogy 
between gravitation and electromagnetism is considered and realized: 
now a practically perfect analogy is attained, without introducing any 
ad hoc hypotheses concerning the gravitational theory. Finally, in chap- 
ter 6 the concluding remarks are made, especially concerning the mod- 
ern development of a canonical formalism in relativistic field theories 
(including the case of gravitation), and studies of the gravity quanti- 
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zation problem. These items are beyond the scope of this paper, but 
they are closely connected with the mathematical methods discussed 
in it, so that we feel it to be worthwhile to mention them. 

1.2 A synopsis of notations of Rieman- 
nian geometry 

A great many of expressions and relations of Riemannian geometry 
will be extensively used below. In existing literature (see for a tra- 
ditional approach, e.g., Eisenhart (1926, 1933, 1972); Schouten and 
Struik (1935, 1938); as to a presentation of general relativity from the 
mathematical viewpoint, see Sachs and Wu (1977)), they are scattered 
in different parts of various publications, and the notations differ sub- 
stantially from one source to another. This is why we give here (usually 
without derivations, though sometimes mentioning original publica- 
tions or principal expositions) quite a few of these geometrical relations. 
The reader may skip this section but look into it merely for clarification 
of some formulae. It is worth emphasizing that we take everywhere the 
speed of light equal to unity (c = 1), the space-time signature being 

(H ); the Greek indices are four-dimensional, running from to 

3, the Einstein summation convention is adopted, including collective 
indices. Symmetrization and antisymmetrization with respect to sets of 
indices are denoted by the standard Bach brackets, (• ■ •) and [• ■ ■] corre- 
spondingly, with the indices in the brackets. Tetrad indices are written 
in separate parentheses, e.g., F^ u y For important details of abstract 
representation of the tensor calculus and its applications to field theory, 
see (Israel 1970, Ryan and Shepley 1975, Eguchi, Gilkey and Hanson 
1980, Choquet-Bruhat, DeWitt-Morette and Dillard-Bleick 1982, von 
Westenholz 1986). 

We admit that the tensor algebra in an arbitrary basis and in ab- 
stract notations, is commonly known. Hence we begin with Cartan's 
forms. To make notations short, we shall sometimes employ collective 
indices, e.g. in the basis of a p-form (cf. Mitskievich and Merkulov 
(1985)): 

dx a := dx ai A • • • A dx ap . 
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Here exterior (wedge) product is supposed to be an antisymmetrization 
of the tensorial product, 

dx ai A • • • A dx a " := dx [ai <g> • • • ® dx a " ] . 

Similarly, we shall write tetrad basis of a p-form as 6^ . Then a p-form 
a can be written as 

a = a a dx a = a (a )# (a) , 

while 

a A = (-l) pq (3 A a (1.2.1) 

where p = deg a, q = deg /3 (degrees of the forms a and (3). 

The axial tensor of Levi-Civita is defined with help of the corre- 
sponding symbol, 

E afh f, = (-g) 1/2 e a ^, E a ™ = -{-g)- 1/2 e aPlP . (1.2.2) 

In our four-dimensional (essentially, even-dimensional) manifold, this 
tensor has the following properties: 

E ag = {-lf^E ga = (-iyE ga , 

E ag E b ° = -p!(4 -p)\5 b a , #a = #b = P = 4-#9- 

While using here collective indices, we denote the number of individual 
indices contained in them, by #. The Kronecker symbol with collective 
indices is totally skew by a definition, 

ft = ■■= <\ • • • € P pV 6 la 6 u] = C (1-2.3) 

so that A a 5 b a = A b (let A a be skew in all individual indices). Dual 
conjugation of a form is denoted by the Hodge star * (see Eguchi, 
Gilkey and Hanson (1980)) before the form. The star acts on the basis, 

*dx a := - — ^——E a a dx 9 , *a := a a * dx a . 
(4-p)! 9 

Thus * * a = (— l) p+1 a (here even-dimensional nature of the manifold 
is essential). It is worth stressing that the sequence of the (lower and 
upper) indices is of great importance too, e.g. it was not the question 
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of a chance when we have written E a g and not E g a (which in fact is 
equal to (-l) p E a g , cf. the formula after (1.2.2)). 

The usual dual conjugation of a bivector will be also denoted by 
a star, but this will be written over (or under) the pair of indices to 
which it is applied, e.g. 

p K x ._ Ie kX ^F^, (1.2.4) 

thus F **p = —F a p. It is self-understood that F K \ is skew (we called it 
bivector for this reason), so that the dual conjugation does not lead to 
any loss of information. Our choice of place for the star is justified here 
by the convenience in writing down the so-called crafty identities where 
the stars are applied to different pairs of indices (see e.g. (Mitskievich 
and Merkulov 1985)): 

n A ;„ = -v„ kX - \s;Xr aT + s;:v aT ^ + s^v^. (1.2.5) 

For special cases and when a contraction is performed, one obtains 
in electrodynamics (cf. (Wheeler 1962, Israel 1970, Mitskievich and 
Merkulov 1985)) 

F/ivF Xu - F; u F": = ~ 5^ F aT F aT , (1.2.6) 



p * jpXv _ \ 6 \p * par_ (1.2.7) 

For the Levi-Civita axial tensor one has 

pnX -pK\ * = _2S"£ (1.2.8) 

for the Riemann-Christoffel curvature tensor, 

? = -Raft* + ~ 4 <^l (1-2-9) 
and for the Weyl conformal curvature tensor, 

C a y s = C a ^. (1.2.10) 
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By the way, a repeated application of the crafty identities to a construc- 
tion quadratic in the curvature tensor (with subsequent contractions, 
so that only two indices remain free ones), leads to the well known 
Lanczos identities, 

R a ^R a ^ u - \R aPl5 R aPlS g^ - 2R^R a(S + R a(3 R^g^- 

-R%R au + R RT - ^flV = (1.2.11) 

(see our derivation of these identities (Mitskievich 1969) which differs 
radically from that by Lanczos (1938) who used less straightforward 
integral relations). Note that the invariant R a p^sR a ^ 7 / is known as the 
Bach-Lanczos invariant, or the generalized Gauss-Bonnet invariant (the 
first of Lagrangians of the Lovelock (1971) series reduces to it); like the 
electromagnetic invariant F a ^F a f , it represents (in a four-dimensional 
world) a pure divergence, so that it does not contribute to the field 
equations in four dimensions. Cf. also DeWitt (1965). 

Scalar multiplication of Cartan's forms is realized with help of con- 
secutive dual conjugations: 

in particular, 

*(dx x A *dx») = -g x ». 
The last relation is equivalent to 

*(E\ p>7 dx x * pa ) = E\ pu E x ^ = -3! g x ». 

The covariant differentiation axioms (see e.g. Ryan and Shepley 
(1975)) are: 

(1) V V T has the same tensor valence as T, v being a vector. 

(2) V„ is a linear operation: 

V V (T + U) = V V T + V V U, 



^au+bvT = aV u T + bV v T. 
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(3) The Leibniz property (distributive property when product expres- 
sions are differentiated) holds, including the contraction operation. 

(4) Action of V on a function (a scalar; the concept of scalar is some- 
what different for the covariant differentiation denoted by a semicolon) 
is V„/ = vf, where v = v a d a . 

(5) Metricity property: V v g = (for a further generalization of the 
geometry, the non- metricity tensor appears on the right-hand side). 

(6) Zero torsion axiom: 

V u v - V v u = [u, v] + T, T = 0. 

(in a generalization of the geometry, the torsion tensor T becomes dif- 
ferent from zero). 

These axioms are most simply realized when the connection coeffi- 
cients are introduced, 

v* (a) *cfl = rjg (a) x (7); v* (a) 6>™ = -rg (a) 0<« 

(only one of these relations is independent). When the structure coef- 
ficients of a basis are introduced with help of 

[X(a),Xffl] = C (a)(/3) X (l)' C (2)(P) = r (J)(a) - r (a)(/3) (1-2.12) 

(for a holonomic, or coordinated basis, = 0; as usually, we write 
in individual parentheses only the tetrad and not coordinated basis 
indices), the general solution satisfying the whole set of axioms, from 
1 to 6, is 

r (a)(/?) = 29 {im [ x {p)9{a){8) +X {a) g mS) -X (j )S( a)M ] 

+ c £m9(«m9 {eKl) + 4ii(fl(flff (()(7) ]- (1-2-13) 

In a coordinated basis (where the basis vectors and covectors are 
usually written as d a and dx a ) the connection coefficients reduce to the 
Christoffel symbols (which are written without parentheses before and 
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after each index, and which are symmetrical in the two lower indices). 
Orthonormal bases and those of Newman-Penrose (see Penrose and 
Rindler (1984a, 1984b)) are special cases of bases whose vectors have 
constant scalar products (i.e., the corresponding tetrad components of 
the metric are constant). For such bases the connection coefficients are 
skew in the upper and the first lower indices; they are sometimes called 
Ricci rotation coefficients. 

Often (in the elder literature, always) covariant differentiation is de- 
noted by a semicolon before the index. Then, according to the covariant 
differentiation axioms (in particular, the axiom 1, however paradoxical 
it would appear), it leads to an increase by unity of the valence of the 
corresponding tensor. As Trautman (1956, 1957) has remarked, such a 
covariant derivative can be defined as 

T a . a :=T a , a + T a \lKa (1-2-14) 

where the coefficients T a \ T a determine the behaviour of the quantity T a 
(which could be not merely a tensor, but also a tensor density or some 
more general object) under infinitesimal transformations of coordinates: 

x '» = x ^ + eC(x), 

5T a := T' a (x') - T a (x) = eT a \lC, T - (1-2-15) 
It is clear that the coefficients T a \ T a possess a property 

(SaTbjYa- = Sa^Tb + SaT^ 

for tensor product S®T, which is closely related to the Leibniz property 
of differentiation. The definition of covariant derivative (1.2.14) may 
be deduced from the concept of Lie derivative, 

£(T a := T a , a C - T a \lC,r = Ta;aC - T a |^. T , (1.2.16) 

which can be written for objects of the most general nature (leaving 
the limits of tensors completely) as 

£(T := e~ 1 (T(x) - T'(x)) (1.2.17) 
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(cf. Yano (1955)). It is important that the both tensors in the paren- 
theses in (1.2.17), have the same argument (non-primed one). This 
leads to a possibility to extend the operation T a . a even to quantities 
of the type of a connection. Without going into details we would note 
that in this sense = R X ( pv ) a . 

It is convenient to define the concept of curvature in the Riemannian 
geometry by introducing the curvature operator (see e.g. (Ryan and 
Shepley 1975)), 

R(u, v) := V U V V - V V V U - V [uM . (1.2.18) 

This operator has the following five properties: 

1) It is skew in u and v. 

2) It annihilates the metric, R(u,v) g = 0. 

3) It annihilates a scalar function, R(u, v) f = 0. 

4) This is a linear operator (cf. the axiom 2 of the covariant differenti- 
ation). 

5) The Leibniz property holds for it (this fact is quite remarkable since 
the curvature operator involves second, not first order differentiation). 

The curvature tensor components (which are scalars from the point 
of view of the V-differentiation) are determined as 

R {K \x )M{ ») := 9^ ■ R(X M ,X {U) )X {X) . (1.2.19) 

These are the coefficients (components) of decomposition of the curva- 
ture tensor with respect of the corresponding basis, and in the general 
case they are 

K (/3)(7)(<5) - X M l (f3)(S) ~ A ( 5 ') 1 (/3)(7) + 1 mS) L M(7)~ 

1 CSXt) 1 (e)(«5) ^(~/)(6) L 03)( e )- Jl.^Uj 

It is worth remembering that some authors use a definition of curvature 
with the opposite sign, e.g. Penrose and Rindler (1984a,b) and Yano 
(1955) (see also a table in (Misner, Thorne and Wheeler 1973)). The 
structure of the curvature operator yields a simple property, 

R\^w x = w K , v ,,-w K ,,, v . (1.2.21) 
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It is also easy to show that simple algebraic identities hold, 

RnXfiv R[k\] [/if] R^lukXi 

as well as 

if K[ A H = = i^V (1.2.22) 

(this last identity bears the name of Ricci). Important role is also 
played by the differential Bianchi identity, 

RkX\p,v£] = = R a f KX .p- (1.2.23) 

The Ricci curvature tensor we define as R^ = R Ufl := R x fiu \, cf. Eisen- 
hart (1926, 1933), Landau and Lifshitz (1971), and Misner, Thorne 
and Wheeler (1973). There are also introduced the scalar curvature, 
R := -R°, and the Weyl conformal curvature tensor, 

C K \^ U := R K \^ + -^(g^Rxu + gxuR^ — gx^Rnu — gnpRxn) 



-qR(9ku9xv - 9kv9x^)- (1.2.24) 

The latter (taken with one contravariant index, and in a coordinated 
basis), does not feel a multiplication of the metric tensor by an arbitrary 
function, so it is conformally invariant. 

In the formalism of Cartan's exterior differential forms, see (Israel 
1970, Ryan and Shepley 1975), the operation of exterior differentiation 
is introduced, 

d:=6^V X(a) A, (1.2.25) 

where A is the discussed above wedge product (the operation of exterior 
multiplication), so that dd = 0. The operation d has its most simple 
representation in a coordinated basis (since the Christoffel symbols are 
symmetric in their two lower indices). It is however especially conve- 
nient to work in an orthonormalized or Newman-Penrose basis when 
the connection 1-forms 

._ pW ail) 
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are skew-symmetric. In the general case they are found as a solution 
of the system of Cartan's first structure equations, 

d0(«) = - w W (ffl AQW, (1.2.26) 

and relations for differential of the tetrad metric, 

dg(a)(l3) = W(a)(/?) +W(/?)(a)- 

These are just 24 + 40 equations for determining all the 64 components 
of the connection coefficients, so that the solution of this set of equations 
should be unique. It is easy to accustom to find the specific solutions 
without performing tedious formal calculations; then the work takes 
surprisingly short time. Since reference frame characteristics (vectors 
of acceleration and rotation, and rate-of-strain tensor) below will turn 
out to be a kind of connection coefficients, this approach allows to 
calculate promptly also these physical characteristics, though for the 
rate-of-strain tensor the computation will be rather specific. 

In their turn, the curvature tensor components are also computed 
easily and promptly, if one applies Cartan's second structure equations, 

(a) (/3) = M<*)w + ^ W (7) A ^ (7 V (1-2.27) 

Here the calculations are much more straightforward than for the con- 
nection 1-forms, and they take less time. It remains only to read out 
expressions for the curvature components according to the definition 

n (a) (fl ^ (Q Wx^ (7) A0 (5) - 

The Ricci and Bianchi identities read in the language of Cartan's forms 
as 

nW M A^ = (1.2.28) 

and 

dQ {a) (f3) = fi (a) {7) A cj (7) (/3) - cu (a) (7) A fi {7) (/3) (1.2.29) 
correspondingly. 

Let us give here a scheme of definition of transports. Denote first 
as S T v :— v\q — v\p^q the change of a vector (or of some other object) 
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under its transport between points P and Q (see (Mitskievich, Yefremov 
and Nesterov 1985)). Then for the parallel transport 

and for the Fermi-Walker transport (cf. Synge (I960)) 
S FW v 



dX 



W u v + r][(v ■ W u u)u - (v ■ u)V u u] =: V F u w v (1.2.31) 



where u a = dx a /d\ is the tangent vector to the transport curve (for 
the Fermi- Walker transport, it is essential to choose canonical param- 
eter along the transport curve in such a way that the norm u • u — r\ 
would be constant and non-zero — so the curve itself should be non- 
null for the Fermi-Walker transport). Using the definitions (1.2.14), 
it is easy to generalize the transport relations to arbitrary tensors or 
tensor densities; so for the Fermi- Walker transport 

%^ = T a , a u a + V T a \l(u T . a u ff - u T u a , a )u a . (1.2.32) 
d\ 

From (1.2.31) a definition of the Fermi- Walker connection seems to 
be directly obtainable, together with the corresponding covariant FW 
differentiation operation which should lead also to the FW-curvature 
{cf. (Mitskievich, Yefremov and Nesterov 1985)): 

r™{ W , V ) : = vrvr - vrvr - v^. (1.2.33) 

The reader may find that in this generalization a difficulty arises: the 
"derivative" V„ w does not possess the linearity property. This trouble 
is however readily removable if we begin with introduction of another 
type of derivative, namely (see below) which possesses all necessary 
properties of a covariant derivative. 

We give here the (principal and equivalent to each other) forms of 
the geodesic equation for a canonical parameter: 

V u u = 0, (1.2.34) 



u A *du = 0, 



(1.2.35) 
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^ + r ?4^-°' (1-2.36) 

, doc 1 doc doc , r\ (~\ i \ 

d\ i9a ^ = 2 9 — ~dX ~dX ■ (L2 - 37) 
A geodesic may be, of course, also null (light-like), but if it is time-like, 
one has d\ = ds. If a non-canonical parameter is used, an extra term 
proportional to u should be added in the equations (1.2.34) - (1.2.37), 
with a coefficient being arbitrary function of the (non-canonical) pa- 
rameter. 

Below we shall make use of the just mentioned generalization of 
the FW-derivative (1.2.31) along a time-like congruence with a tangent 
vector field r normalized by unity (the monad). It is convenient to 
denote such a derivative as V^, 

W U T = V U T + TX {r^ - T^, a ) u a B a , (1.2.38) 

where B a is the (coordinated) basis of T: T = T a B a . Then V^ w = V^. 
The definition (1.2.38) can be also written as 

T a; M = T a ,, + 7X(r% ;M - r%r p ). (1.2.38a) 

Consider now the properties of this r-derivative when it is applied to a 
scalar function /, the vector r, the metric tensor g, and the construction 
b = g—T®r (the three-metric or the projector in the monad formalism): 

V£/ = «/, V^r = V^ = V^6 = (1.2.39) 

(action of a r-derivative on a tensor is understood in the sense of 
(1.2.38); u is considered as a partial differentiation operator in uf). The 
corresponding r-curvature operator is denoted as R T (w, v) (cf. (2.4.10); 
see also (Mitskievich, Yefremov and Nesterov 1985)). 

Returning to the Lie derivative, it is worth mentioning the relations 

%^ = ^ + &^ (1-2-40) 



- i x -,w - CR x ^k = £ A ( ; ^) - CR X ( P u)k- (1.2.41) 
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A vector field satisfying the equation 

Z„.r I 0. (1.2.42) 

is called the Killing field. Its existence imposes certain limitations on 
the geometry, — existence of an isometry: £^g (liV = (see (Yano 1955, 
Eisenhart 1933, Ryan and Shepley 1975)). The Lie derivative plays 
an important role in the monad formalism, as we shall see below; this 
is why a version of this formalism well adapted to description of the 
canonical structure of the field theory, is called the "Lie-monad formal- 
ism" (see (Mitskievich and Nesterov 1981, Mitskievich, Yefremov and 
Nesterov 1985)). 

Alongside with exterior differentiation d which increases the degree 
of a form by unity, in Cartan's formalism there is employed also another 
differential operation which diminishes this degree by unity, 

5:=-*d*. (1.2.43) 

In fact this is the divergence operation, 

5(a^ h dx^) = pa^dx\ p = #h + l, (1.2.44) 

which does not however fulfil the Leibniz property: when applied to a 
product, e.g., of a function and a 1-form, fa, it results in 

5(fa) = f5a — *(df A *a), 

while it acts on a function (0-form) simply annihilating it, 5f = 0. One 
might call this a generalized Leibniz property. Of course, similarly to 
the identity dd = 0, there holds also the identity 55 = 0. A combination 
of the both operations gives the de Rham operator (de-Rhamian), 

A := d5 + 5d = (d + 5) 2 (1.2.45) 

(remember divergence of a gradient in Euclidean geometry). A form 
annihilated by the de-Rhamian, is called harmonic form (remember 
harmonic functions and the Laplace operator). In contrast to 5, the 
de-Rhamian acts non-trivially on 0-forms, and in contrast to the ex- 
terior differential d, it acts non-trivially on 4-forms too (also in the 
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four-dimensional world). There are introduced also the following clas- 
sification elements of forms: an exact form is that which is an exterior 
derivative of another form; a closed form is that whose exterior differ- 
ential vanishes (an exact form is always closed, though the converse 
is in general not true); a co-exact form is a result of action by 5 on 
some form of a higher (by unity) degree; a co-closed form is itself iden- 
tically annihilated by 5. In (topologically) simplest cases it is possible 
to split an arbitrary form into sum of a harmonic, exact, and co-exact 
forms (see Eguchi, Gilkey and Hanson (1980), Mitskievich, Yefremov 
and Nesterov (1985)). 

For an arbitrary form a, 

A« = a a .^dx a - pR\^a^ h \ T a dx vh , #a = p. (1.2.46) 

The number of individual indices entering the collective index a of the 
components of the form a, may be any from to 4, while the quantity 
a a \^ is determined by (1.2.15). The expression (1.2.46) can be reduced 
identically to 

Aa = (pO-^cvm* +pR°a auh + V ^a^daT* (1.2.47) 

where a = fj,uh, j^a = 2 + #h = p, and p may be equal also to zero. 

We compare here explicitly the results of action of de-Rhamiam on 
forms of all degrees in a four- dimensional world: 

/: A/ /,/". (1.2.48) 
A = A a dx a : AA = {A a .» + R a a A a )dx a , (1.2.49) 



F = l -F aP dx a(3 : AF = ^{F a p,^+2RlF aP +R^ aP F aT )dx a ^ (1.2.50) 



^T al3y dx a ^ : AT = ^(T af3r ^ + SR a a T af3l + 3R° T a(3 T aT ^)dx a ^ 
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(A^ + R°A a ) * dx a , 



(1.2.51) 



V = -V aPl5 dx a ^ 5 : AV = -V a ^dx a ^ s = f$ * 1. (1.2.52) 



From the last two expressions it is clear that there exists a symme- 
try with respect to an exchange of a basis form to its dual conjugate 
with a simultaneous dual conjugation of components of the form (com- 
pare (1.2.48) and (1.2.52), (1.2.49) and (1.2.51)). Here we made use of 
obvious definitions 




T al3y E al3 "' s , / — ——V al3y sE' 



alongside with the easily checkable identity, 



{2R a a E a ^i, + 3R <TT a pE aT1 s)E a 7 = 0. 



Chapter 2 



Reference frames' calculus 



2.1 The monad formalism and its place 
in the description of reference frames 
in relativistic physics 

It would be completely erroneous to automatically identify the concepts 
of a system of coordinates and of a reference frame. Such an identifi- 
cation is often characteristic for those physicists who are accustomed 
to work in Cartesian coordinates and who, moreover, try to reduce (for 
sake of "simplicity") problems — mainly, educational ones — to two- 
dimensional diagrams dealing with one temporal and another spatial, 
coordinates. But when more than one spatial dimension (plus extra 
one, temporal dimension) are considered, a question arises whether it 
is possible to relate to different reference frames systems of coordinates 
which differ one from the other by a purely spatial rotation (in the 
simplest case when the both systems are Cartesian). The more such 
a problem becomes obvious when we pass from a Cartesian to polar 
system of coordinates. 

Let us call to mind the Newtonian mechanics (see, e.g., Newton 
(1962), Khaikin (1947), H. Goldstein (1965), Schmutzer (1989), as well 
as Schmutzer and Schiitz (1983)), and think it over, in which cases there 
appear, for example, the forces of inertia (we do not intend to take part 
in scholastic discussions over the reality of those). What is clear is that 
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such forces are originated by the state of motion of the reference frame. 
But if we consider two reference frames at rest, one with respect to the 
other, — could any difference between them be physically perceivable in 
course of some experiments, at least in principle? If this would be the 
case, it would be possible to experimentally determine, what kind of 
coordinate frame (e.g., Cartesian or polar one) is used by the observer, 
or, at least, which of the relatively resting bases corresponds to the 
latter. It is however clear that experiments are performed irrespective 
to the choice of purely spatial coordinates (the corresponding system 
of coordinates is chosen only at the stage of mathematical description 
of physical effects, and the experimental devices — even if they are 
used for measurement of lengths or angles — have no direct relation 
to this choice). The same relates also to the choice of a purely spatial 
basis. Thus transcriptions of all expressions and quantities supposed to 
describe physical effects, have to be performed invariantly with respect 
to the choice of spatial coordinates or spatial basis. This demand is 
not as mild as it may seem to be at the first sight, and it yields a set of 
geometrical identities to be imposed on the corresponding congruences. 

Hence we see that a transition between different reference frames 
should basically reflect their relative motion. The simplest special case 
of such a motion is that which relates two non-identical inertial refer- 
ence frames both in Newtonian and relativistic mechanics. 

It is also obviously possible to an observer or a system of observers 
to move together with a frame of reference, locally or in a space-time 
region (globally). When an object moves together with a reference 
frame, it is geometrically identified with the latter, so that if this is 
a spatially extended object, the world lines of its mass points form a 
congruence. Such a congruence presents a complete characterization of 
the reference frame, and with any reference frame a conceptual object 
(thus a test one) of the above type can be associated which is called the 
body of reference. Since it models a set of observers and/or their mea- 
suring devices, but not photons, the reference frame congruence should 
be time-like. A transition to another reference frame means merely a 
change of such a congruence, with the corresponding and automatic 
change of the test reference body. In order to avoid mathematical dif- 
ficulties, one has to necessarily use congruences, since otherwise there 
would appear a danger of either mutual intersections of the world lines 
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(which mathematically manifests itself as a non-physical singularity), 
or in the space-time there would arise bald spots whose boundaries will 
represent singularities too. 

The simplest way to describe a reference frame is to identify the 
congruence of the corresponding reference body mass-points with a con- 
gruence of the time coordinate lines. The latters should be, of course, 
time-like, which is not always the case for the so-called time coordinate 
(e.g., for the Schwarzschild space-time in the curvature coordinates, the 
t coordinate lines become null on the horizon and space-like inside it), 
thus we come to a restriction concerning the choice of coordinates ad- 
equate for the reference frame description of physical situations. Here 
a strange (from the point of view of naive psychology) fact arises: it is 
easy to see that any transformation of the time coordinate with fixed all 
other (spatial) coordinates, can at most change the lines of the latters, 
but not the time coordinate lines, i.e. we remain with the former refer- 
ence frame if determined through the time coordinate description. Such 
transformations are called the chronometric ones, and the invariance of 
a reference frame with respect to them gives a basis for Zel'manov's 
formalism of chronometric invariants (Zel'manov 1956, 1959) being of 
importance primarily for cosmology (see also (Mitskievich and Zaharow 
1970, Mitskievich 1969)). 

Another formalism, that of kinemetric invariants, was also proposed 
and developed by Zel'manov (1973). It is based on a description of ref- 
erence frames with help of families of space-like hypersurfaces in the 
space-time. We confine ourselves here to the corresponding reference 
since, as also the formalism of chronometric invariants, this is in fact 
a special choice of the system of coordinates in the monad formalism 
extensively discussed below. This latter was primarily worked out by 
Zel'manov (1976); see also good reviews of all these three formalisms 
in Zel'manov and Agakov (1989), Massa (1974a,b,c), and partially in 
Shteingrad (1974), Mitskievich (1972), Ivanitskaya, Mitskievich and 
Vladimirov (1985, 1986). It should be noted that the formalism of 
chronometric invariants, with some modifications (which bring it near 
to the monad formalism), was independently developed by Cattaneo 
(1958, 1959, 1961, 1962) and Schmutzer (1968), while some applications 
to the study of physical effects are given in Schmutzer and Plebanski 
(1977) and Dehnen (1962). In its turn, the formalism of kinemetric 
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invariants represents a refinement of the earlier formalism of Arnowitt, 
Deser and Misner (ADM) (see Misner, Thorne and Wheeler (1973), as 
well as Mitskievich and Nesterov (1981), Mitskievich, Yefremov and 
Nesterov (1985)). This ADM formalism proved to be of importance 
for realizing the canonical formalism and the quantization procedure 
in gravitation theory. For applications of an approach similar to the 
monad formalism, to astrophysical effects, see Thorne and Macdonald 
(1982), Macdonald and Thorne (1982), Thorne, Price and Macdon- 
ald (1986). An extensive review of methods of description of reference 
frames see in a monograph by Vladimirov (1982) who himself has made 
a substantial contribution to this field, and also in a recent paper by 
Jantzen, Carini and Bini (1992). 

During many years the principal role in description of reference 
frames in general relativity, was played however by the tetrad formal- 
ism which influenced much the very style of thinking. In the case of 
an orthonormalized tetrad, we are dealing not only with a time-like 
congruence which corresponds to four of the sixteen components of the 
tetrad (the field of tangent vector to the congruence), but with three 
(spatial) congruences more. One does not usually care much about the 
invariance of expressions which have a physical meaning, although this 
obviously contradicts the conclusions drawn above. If we take such a 
care, we come to the monad formalism with auxiliary use of an arbi- 
trary triad in the subspace orthogonal to the monad which does not 
practically change our further conclusions. The tetrad formalism was 
considered by many authors; for an orientation in the corresponding lit- 
erature, see Ivanitskaya (1969), Mitskievich (1969) and, in connection 
with gravitational effects, Ivanitskaya (1979). 

Tetrad bases are not always directly applicable to description of ref- 
erence frames, e.g. in the case of a Newman-Penrose basis (see Penrose 
and Rindler (1984a, 1984b) where also the use of spinors in constructing 
bases is considered). In some calculations, it proved effective to combine 
together the Dirac matrices and Cartan forms when gravitational fields 
and reference frames are studied (Mitskievich 1975). General methods 
of the tetrad formalism function well in the coordinates-free approach 
to tensor calculus, in particular in the Cartan forms formalism (see 
section 1.2). 

We came thus to a conclusion that it is namely the monad formalism 
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which is most closely connected with description of frames of reference, 
and this formalism is presented in the next Chapter. A question is 
however still worthwile: is it really always indispensable to use explicitly 
some definition of a specific reference frame when observable physical 
effects are calculated from the point of view of a non-inertial observer? 
And in what classes of problems of theoretical physics, and at which 
step of calculation procedure, the formalisms of description of reference 
frames may be of practical use? 

When local effects are studied, this means that the observer may 
be idealized as a single point-like test mass while he receives all outside 
information through signals propagating with the fundamental velocity 
(electromagnetic and other signals), as well as detecting all other phys- 
ical and geometric factors which are perceivable on his world line only. 
Then this world line will be the only essential representative of the 
state of motion of the observer in the calculations, so that any exten- 
sion of the reference frame to the outside will be absolutely unnecessary. 
Moreover, if the magnitude of some predicted physical effect would to 
any extent depend on the concrete way of this extension, this effect 
could not be considered as a correctly determined physical one, since 
there should be an invariance of all observables (for our local observer) 
with respect to an arbitrary choice of the reference frame outside the 
observer's world line, if only this line would, with a proper smooth- 
ness, belong to the congruence which globally describes the pertinent 
reference frame. It is obvious that there exists infinite and continuous 
multitude of such reference frames mutually coinciding on a single world 
line. An example of a calculation of a physical effect when no extension 
of the reference frame outside a single world line of the observer should 
be considered, is the case studied in section 3.4, a united description 
of the gravitational red shift effect and the Doppler effect. These two 
effects are inseparable when the space-time is non-stationary, and their 
distinction is a mere convention. This conventionality manifests itself 
in the mentioned case in the following two interpretations of the effect. 
When the both time-like world lines (those of the emitter of radiation 
and of its detector) are considered to belong to one and the same con- 
gruence of a reference frame, the both emitter and detector are at rest 
by definition, and the effect has to be interpreted as a purely gravita- 
tional one. The distance between the two objects is however continually 
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changing; in a similar situation in cosmology, the non-stationarity of 
space-time completely prevents domination of any of the alternative 
interpretations! Another approach consists of inclusion in the reference 
frame congruence of the observer's (detector's) world line only; then 
the radiation source will be in motion with respect to the frame, so 
that the effect would be (partly) due to the Doppler shift. An interpre- 
tation of this effect as a completely Doppler one, is also possible since 
the space-time is stationary; the reference frame congruence is then the 
Killing one, although it fails to be time-like in the ergosphere. But 
outside the latter, one has to consider the effect as being due to mo- 
tions of the both emitter and detector of radiation. As to the time-like 
Killing congruence, if it is rotating, one usually connects it with the 
dragging phenomenon. Such an objectively determined reference frame 
is discussed in section 3.5. 

There exists however a wide class of still other problems in which 
employment of a concrete global reference frame is perfectly accept- 
able and desirable, while the predicted in this way effects do depend 
on its specific choice. These are the problems in which the interacting 
physical factors possess spatial extension, e.g. the case of charged per- 
fect fluid moving in electromagnetic and gravitational fields. In section 
4.4 we consider such an exact self-consistent problem which moreover 
yields a strange (at first sight) result that for the state of motion under 
consideration, being consistent with behaviour of fields, the non-test 
charged perfect fluid produces in its co-moving reference frame only a 
magnetic, but no electric field whatsoever. Hence its particles, being 
undoubtedly electrically charged, do not mutually interact electromag- 
netically. We show there also that there is no paradox, but this is 
merely an encounter with unusual (for the prevailing inertial mode of 
thinking) specific feature of electromagnetic field in non-inertial, in par- 
ticular, rotating reference frames. We shall discuss there a similarity 
of this situation with that in classical mechanics with its centrifugal 
and Coriolis forces. The concept of a reference frame and employment 
of formalisms of its description is certainly fruitful in the relativistic 
cosmology where the continuous medium filling the universe, realizes 
a privileged (really existing) reference frame. The fact of a non-test 
nature of this medium does not mean that the reference frame per- 
turbs in this case geometry, this is merely a specific — co-moving - 
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frame, which is objectively singled out of all other frames by the very 
statement of the problem. 

We do not (and cannot) give here an unambiguous and precise an- 
swer to the question what problems demand application of the concept 
of reference frame, and for what problems this is unnecessary. Instead 
we confine ourselves to certain examples, since it would be a scholastic 
task, to foresee every special case of application of reference frames in 
the future. As to introduction of this concept at an initial stage of every 
calculation, this clearly cannot be practical, since symmetry of a prob- 
lem may more easily be connected with a system of coordinates and 
not with a reference frame, so that the adequate choice of coordinates 
is more crucial for the most effective treatment of problems (e.g., when 
Einstein-Maxwell fields are derived). When however one turns to in- 
terpretation of the obtained solutions, introduction of reference frames 
becomes quite opportune, and it may play a decisive role. These consid- 
erations bring us to a better understanding of the fact that the first of all 
complete and practical methods of description of reference frames was 
Zel'manov's formalism of chronometric invariants (the quotient system 
of coordinates co- moving with a frame of reference), and this was pri- 
marily done in problems of relativistic cosmology. In the same context, 
certain problems connected with the Noether theorem were also con- 
sidered by Mitskievich and Mujica (1967) and Mitskievich and Ribeiro 
Teodoro (1969). 

The concept of reference frame should be probably essential also at 
starting points of development of new theoretical ideas. As an example 
the idea of zero-point radiation may be mentioned which seems to be 
still at a relatively early stage of its development. We only recall that 
the zero-point radiation possesses such a spectrum that under Lorentz 
transformation, it retains its form (while a black-body radiation would 
alter its temperature), but a transition to a uniformly accelerated frame 
yields in addition to the former zero-point radiation, also a black-body 
thermal one (Boyer 1980). The reader should not be perplexed by the 
fact that the integral energy of the zero-point radiation diverges (it is 
usually treated in the spirit of renormalization ideas). These items are 
closely related to the concept of the Rindler vacuum, and they accord 
with the quantum theory of black holes and the black hole evaporation 
(Thorne, Price and Macdonald 1986). In quantum physics the concept 
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of reference frame has however not yet found its adequate realization. 
At the same time, its possible application areas are quite far-reaching: 
this is the theory of quantum fields itself (Gorbatsievich 1985) as well 
as motivation and fundamentals of the field quantization (the ADM for- 
malism to be once more mentioned (Misner, Thorne and Wheeler 1973), 
alongside of a new development in this area — Ashtekar's canonical for- 
malism (Ashtekar 1988), some application of which will be sketched out 
in the concluding chapter of this paper). 

As it was motivated above, we have chosen the monad approach 
for description of frames of reference, while monad means a fixed unit 
time-like vector field which is naturally interpreted as a field of four- 
velocities of local observers equipped with all necessary measurement 
devices. In the space-time region where such a reference frame is re- 
alised, we may therefore consider a congruence of integral curves of 
this vector field, i. e. of the world-lines of particles forming a reference 
body. Every local observer represents such a test particle. The classical 
(non-quantum) physics approach admits that acts of measurement do 
not influence the objects and processes which are subjected to the mea- 
surement, in full conformity with our treatment of local observers as 
test particles, — in other words, in conformity with the identification 
of a reference frame and a time-like congruence which, of course, should 
be geometrically admissible in the space-time region under considera- 
tion. Then a transition to another reference frame should be equivalent 
to a new choice of such an admissible congruence, without any change 
of the four-geometry of the world itself. 

Rodichev (1972, 1974) has undertaken an attempt to overpass the 
limits of this approach to reference frames, but we consider his attempt 
as a premature and technically wrong one, though it is worth mak- 
ing an acquaintance with his concepts which were expressed with much 
inspiration. These concepts might become useful for far-reaching gener- 
alizations in the future, such as those connected with quantum gravity. 
As a matter of fact, in quantum theory it is useless to try to introduce 
test particles in the classical meaning, and moreover, the local observers 
cannot be considered there without perturbations of the four-geometry, 
as well as of all involved physical objects. Thus one has to consider a 
radical revision of the concept of reference frame when passing to quan- 
tum physics. Some ideas of such a generalization were outlined by Fock 
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(1971) (see also (Schmutzer 1975)), and some (rather obscure) remarks 
were made by Brillouin (1970), but probably the most powerful insight 
can be found in the early paper by Bohr and Rosenfeld (1933) (see also 
Regge (1958) and Anderson (1958a,b; 1959)). We intend here, however, 
to confine ourselves to classical relativistic physics. 



2.2 Reference frames algebra 

The first steps in the development of the monad formalism were made 
already in the forties and fifties by Eckart (1940) and Lief (1951) on the 
basis of hydrodynamical analogies. In our general relativistic approach, 
this means that in the region under consideration, a time-like four- 
vector monad field r is introduced in addition to the four-dimensional 
metric g, with a normalization condition r • r = +1. Then a symmetric 
tensor, 

b = g-r®r, (2.2.1) 

will be automatically the projector onto the local subspace orthogonal 
to r: 

6(r,-) = 0, b x x = 3, I 

det^ = 0, 6^ = 6J. J [Z - Z - Z) 

As we shall see now, the tensor b plays simultaneously the role of a 
metric in this subspace (the three-dimensional physical space of the 
reference frame). The three-space is then non-holonomic in general 
(namely, when the r-congruence is rotating). 

Consider now an arbitrary four-(co)vector q; its projection onto the 
monad is a scalar, 

(t) 

Q:=q-T, (2.2.3) 

and the projection onto the three-space of the reference frame, a four- 
dimensional (co)vector, 

( q:=b(q,-), (2.2.4) 

which is by a definition orthogonal to the monad and not changing by 
a repeated projection. Hence, 



(r) (3) 

q = Q T+ q 



(2.2.5) 
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Applying the same procedures to another (co)vector, p, one can write 
now 

(t)(t) (3) (3) (t)(t) (3) (3) 

p • q = g(p, q) = p q + p ■ q=:p q - p • q , (2.2.6) 

where we have introduced the notation • for a three-dimensional scalar 
product relative to the reference frame r. The sign before this three- 
scalar product is chosen according to signature of the tensor b, which 
can be written symbolically as 

sign b =(0, -,-,-), (2.2.7) 

it guarantees non-negativity of such a three-scalar square of any vector. 
Now observe that 



p. q = *[( r Ap) A *(r A q)] (2.2.8) 

where the (co)vectors p and q are considered to be arbitrary (four-di- 
mensional), and their projection onto the physical three-space of the 
reference frame is provided automatically by the projecting properties 
of the three-metric. 

Let us apply this simple result to the squared interval using nota- 

(3) 

tions dt and d x for the elements of physical time and three-spatial 
displacement (relative to the reference frame under consideration) cor- 
respondingly according to the definitions of projections, (2.2.3) and 
(2.2.4). These elements are, of course, not total differentials (i.e., they 
are non-holonomic). Thus 

ds 2 = dt 2 - dl 2 , dl 2 := d { x -d ( x= -b(dx, dx), (2.2.9) 

and for propagation of a light pulse (signal), 

ds 2 = (1 - v 2 )dt 2 = v = l, (2.2.10) 

where v = dl/dt is the absolute value of the three- velocity. This result 
is a standard assertion of the universality of the speed of light (but now 
in an arbitrary reference frame and in any gravitational fields). This 
conclusion does not contradict the experimental facts of the time-delay 
of electromagnetic signals when propagated near large masses (e.g., that 
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of the Sun), since in such experiments, the distances were calculated in 
the same sense as in our approach, but the physical time intervals were 
measured by a clock on the Earth, while in the expression (2.2.10), 
the local clocks are supposed to be used which should be localized 
continuously along the path of the signal (this means, of course, merely 
the corresponding re-calculation, as it was the case for the distances!). 

It is clear that the four- velocity of any object can be easily expressed 
through the monad vector and the three- velocity of this object, 

while 

^= ( M ) =(l-t; 2 )- 1 /2. (2.2.12) 

Multy-index quantities (components of tensors) are to be projected 
in each of their indices onto either the monad, or the three-space or- 
thogonal to the latter, so that the whole family of objects of various 
tensor ranks arises (some examples can be found below). 

Let us introduce the three-spatial (and hence three-index) axial 
Levi-Civita tensor, 

e Xlxv := t k E kX ^. (2.2.13) 

Practically, we shall not use this objects, but its existence is obvious, 
and considering it helps to understand more profoundly some expres- 
sions below. Similarly to (1.2.3) we introduce also the skew tensor 
with antisymmetrization in all individual indices forming each of the 
collective ones. Then 

= -0!3!&$; e a ^ = -1!2!6#; \ . 
e aP ve a *> = -2\UbZ; e af3l e a ^ = -3!0!. J 1 ' ' ] 

As it was the case in the flat three-dimensional world for the Levi- 
Civita symbol, the axial tensor (2.2.13) makes it possible to relate rank 
two skew tensors to the corresponding vectors which also belong to the 
three-space of the reference frame. This fact becomes, as usually, the 
basis for definition of the vector product (and, later, the curl operation, 
(2.4.1)). Its symbolic expression reads 

p x q := *(p A r A q). (2.2.15) 
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As in the scalar product (2.2.8), one may take here simply the four- 
vectors p and q without specially projecting them onto the three-space 
of the reference frame, while the resulting vector product will depend on 
the properly projected vectors, and it belongs itself to this three-space. 
The corresponding expression for a mixed triple product is 

n • (p x q) = *(r A n Ap A q). (2.2.16) 

It is convenient to introduce 1-form quantities 

&W=0(«)- r («) r , b (a) = X (a) - T (a) r (2.2.17) 

(the coordinated basis could be also used). Now the three- metric itself 
arises as a three-scalar product, 

*(6< a > A = *(& (/3) A *& (a) ) = 6 (Q) • b i/3) = -b^W (2.2.18) 

while we have, naturally, 

*(M Q ) A *# (/3) ) = A *6^) = -b {a ^\ (2.2.19) 

Other important formulae involving r and b^ are 

dx a/3 = T a T Ab f3 _ T P T /\ b a -\- b a A b 13 (2.2.20) 

and 

dx a<i = T a r A b p - t?t Ab a - e a/3 7 * (r A b r ). (2.2.21) 

Transition from the first to the second relation means that there holds 
identity 

b a A IJ 3 = -e Q/3 7 * (r A 6 7 ) (2.2.22) 

which can be easily checked. 

The monad formalism was first introduced by Zel'manov (1976). 
This monad description of reference frames is connected with the for- 
malism of chronometric invariants (which belongs also to Zel'manov 
(1956)) by a special choice of the system of coordinates. It is impor- 
tant that the both approaches do not involve any restrictions upon the 
choice of a reference frame, except for the standard limitations due to 
geometry and physics in general, but the monad formalism allows in 
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addition also an arbitrary choice of systems of coordinates which are 
then by no means connected to the reference frame. 

Zel'manov's formalism of chronometric invariants is in fact the mo- 
nad formalism written in coordinates co-moving with the monad, i.e. 
such that r M = ^o/v^oO' This means that in this case the monad 
congruence coincides with the time coordinate lines. This coincidence 
does not change under arbitrary transformations of the time coordinate 
(so that this may arbitrarily mix with the spatial coordinates), but the 
new spatial coordinates should be (any) functions of the old spatial 
coordinates only. This special choice of r substituted into all algebraic 
and differential relations of the monad formalism, yields the complete 
structure of Zel'manov's formalism of chronometric invariants. The 
latter may seem to be intuitively more manageable than the monad 
formalism is, but in the reality relations involving chronometric inva- 
riants are more complicated; moreover, in the monad formalism one 
may attain additional simplifications by choosing special systems of 
coordinates which are by no means constrained by an already fixed 
reference frame. This is why we shall not consider below the formalism 
of chronometric invariants in more detail. 

2.3 Geometry of congruences. Acceler- 
ation, rotation, expansion and shear 
of a reference frame 

As we have seen above, when we discussed in section 2.1 the setting of 
reference frames formulation, the crucial role in their description has 
to be played by a study of congruences of the world lines of particles 
forming bodies of reference, i.e. the physical time lines for the chosen 
reference frame. The congruence concept is essential because for the 
sake of regularity of the mathematical description of the frame, these 
lines have not to mutually intersect, and they must cover completely 
the space-time region under consideration, so that at every world point 
one has to find one and only one line passing through it. Exactly the 
same approach is used for description of a regular continuous media, 
i.e. in hydrodynamics of a compressible fluid. Therefore it is worth 
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using here the standard formalism of hydrodynamics (relativistic: Lich- 
nerowicz (1955), Taub (1978), non-relativistic: S. Goldstein (1976)). In 
a line with these ideas we introduce below the concepts of acceleration, 
rotation and rate-of-strain for a frame of reference. 

Differential operations on a manifold can (and for our purposes have 
to) be projected onto the physical time direction of a frame of refer- 
ence and its local three-space with help of the projectors r and b, as 
it was the case for tensor quantities, in particular, exterior forms. It 
is however important to keep in mind that when quantities lying in 
the three-space of a reference frame are being differentiated, new terms 
arise which are non-orthogonal to the monad field. Hence in order to 
obtain after such a differentiation purely three-spatial quantities, one 
has to perform projections not only in the differentiation indices, but 
also in those which stem from the very quantities being differentiated: 
in general, neither the three-projector, nor the monad vector are con- 
stant relative to the differentiation operations, if the latters are not 
reformulated specially for this end. We mean here the operations with 
the properties of the type of (1.2.39), which are to be used for this new 
differentiation. Another operation which is evidently destined for tak- 
ing a derivative with respect to the physical time of a reference frame, 
is the Lie derivative with respect to r; but it does not commute either 
with the covariant components of r, or with the three-projector (both 
covariant or contravariant). Moreover, the Lie derivative, to a greater 
extent than the nabla operator, depends on the choice of norm of the 
vector with respect to which the differentiation is performed (a scalar 
cannot be brought from this vector since it enters the Lie derivative 
under the differentiation sign, see (1.2.16). We have made use of this 
"deficiency" when discussing hypersurface geometry in (Mitskievich, 
Yefremov and Nesterov 1985) — the integrating factor method works 
here in the case of normal congruences only. 

The formalism we treat below, is well fitted not only for the reference 
frames description, but also, to a full extent, in studies of relativistic 
hydrodynamics. In an analogy with four-dimensional hydrodynamics, 
we may write the covariant gradient of monad vector as 

T li[V = T v G li + A vli + D vli , (2.3.1) 

where G is the acceleration (co)vector of reference frame, A the angular 
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velocity tensor or rotation of the frame, and D the rate- of -strain tensor 
(the sum A + D represents the purely spatial part of the gradient, which 
is further separated into symmetric and skew parts, D and A). Let us 
consider now these quantities more closely. 

Acceleration is the most obvious of these constructions. It is clear 
that it should be equal to V T r; a substitution of the expression (2.3.1) 
brings this to 

V r r = tV^ a = G. (2.3.2) 

Since the norm of r is equal to unity by a definition, t^t m;j , = 0, and 
all other terms lie completely in the three-space of the reference frame. 
From the definition of Lie derivative and (2.3.1), we know that 

£ri~n = t^ v t v + t v .^t v = (2.3.3) 

which yields 

£ T b tlv = £ T g jXV - t 11 G v - t u G^ = 2D^ V (2.3.4) 

(see also (2.3.1)). This conforms to the understanding of the rate-of- 
strain tensor as a measure of evolving the spatial scales, these being 
determined by the three-metric b. Finally, rotation of the r-congruence 
is determined as a covector 

u :=^*(r Adr) = *(r AA), A := ^A^dx^ . (2.3.5) 

(from here one-to-one correspondence is obvious between the bivector A 
and axial (co)vector u, the both bearing one and the same information 
about rotation of the reference frame). Hence 

A = - * (r Aw), G = - * (r A *dr). (2.3.6) 

The rate-of-strain tensor (2.3.4) splits naturally into its trace (the scalar 
of expansion, or dilatation) 

© := \r a -, a = ee (2.3.7) 
and the traceless part (the shear tensor), 

<V := D,,,, - -Gb^, (2.3.8) 
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whose square 

a 2 = o^ow ee a^D^ = D, U D^ - ^0 2 (2.3.9) 

invariantly characterizes the presence of shear by evolution of the three- 
space of a reference frame. These concepts were borrowed by the ref- 
erence frames theory from hydrodynamics where they play an impor- 
tant role. They all are important also in theory of null (light-like) 
congruences often used in the classification of gravitational fields by 
principal null directions (the Petrov types) and generation of exact 
Einstein-Maxwell solutions (see Kramer, Stephani, MacCallum and 
Herlt (1980)). When a null complex (Newman-Penrose) tetrad is con- 
sidered, some of the above relations change their form. 

Consider now exterior differential of the monad covector, 

dr = rAG + 2A (2.3.10) 

which can be rewritten as 

dr = -G At + 2*(u At) (2.3.11) 

(cf. (3.1.16) for the electromagnetic field case). This relation is a part 
of Cartan's first structure equations adapted to a monad basis which 
can be supplemented (without leaving the coordinated basis approach) 
by quantities similar to (2.2.17), a covector 

ff = dx» - t^t, (2.3.12a) 

and vector 

b„ = dp - t^t. (2.3.12a) 

One has to keep in mind that the free index \i appearing in these 
expressions, does not contribute to the tensor properties of the involved 
quantities from the point of view of the coordinate-free methods (the 
operations V and exterior differentiation), but it spoils the properties 
of the components of these quantities so that Christoffel symbols enter 
such expressions, e.g. 

dff = -2rM + [r"G„ - (A/ 4 + AS) + r^r A ]6" A r (2.3.13) 
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Here these symbols play role of scalar functions (from the point of view 
of exterior forms formalism). One can however get rid of them by 
introduction of a differential operation 

(b) 

d := b a X7 b A = b a X7 da A = -rV T A . (2.3.14) 
Making use of the relation 

r AV/ = K A r¥ + r"G)Ar, (2.3.15) 

we obtain finally 

{ d W = -2r"A - (A/ + A/)b v A r (2.3.16) 

Here (cf. (2.3.5)) it is possible to write A = (1/2) A^ A 6". The price 
one has to pay for this simplification, is incompleteness of the system 
of structure equations (2.3.10) and (2.3.16), which however does not 
prevent obtaining from them complete information about the quantities 
G, D and A. Let us simultaneously write down relations analogous to 
(2.3.2): 

V b ^ = (D^ + A^)b v , ( 2 - 3 - 17 ) 
V r 6 M = —TpG - G M r + r^r^ A , (2.3.18) 
V 6a A = (^A " r v T$b K - r WiA r)6j. (2.3.19) 



2.4 Differential operations and identities 
of the monad formalism 

When a reference frame is not taken into account and only the four-di- 
mensional nature of the universe matters, the role of curl is played by 
exterior derivative, and role of divergence by the operation 5 (1.2.43). 
With help of these operations the four- dimensional de-Rhamian is then 
built which represents a generalization of the Laplacian (d'Alembertian) 
operator to the theory of differential forms. In order to introduce anal- 
ogous concepts in the scope of monad description of reference frames, 
it is sufficient to recall their definitions in three-dimensional Euclidean 
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geometry where, e.g. curl is considered as a vector product in which 
the first factor is the nabla operator. Here we shall base on the vec- 
tor product expression (2.2.15), but the monad vector r will be not 
influenced by the exterior differentiation, 

curl ( a = *(r A d ( a ). (2.4.1) 

In this definition, curl is acting on a covector which is already lying in 
the three-space of the reference frame. Remember that in the defini- 
tion (2.2.15) it was unnecessary to consider any of the covectors taking 
part in the vector multiplication, to be projected onto the three-space: 
they were projected automatically by virtue of (2.2.15). But in the 
case of curl this kind of definition leads to inclusion of an additional 
term compensating the part proportional to r of the covector being 
differentiated. It is easy to see that then 

(t) 

curl a = *(r A da) - 2 a cu. (2.4.2) 

It might seem that in this expression the overall sign should be chosen 
in the opposite sense (cf. Mitskievich, Yefremov and Nesterov (1985)), 
but the expression (2.4.2) is perfectly correct: as a matter of fact, 
the Euclidean expression for curl being a contravariant vector, it con- 
tains the differentiation operation (nabla) which is treated there as if 
it yielded directly a contravariant index. It is however the simplest to 
compare the curl components (2.4.1) brought to the contravariant level, 
and its components in the standard Euclidean vector calculus. 

Let us turn now to the three-divergence operation. The natural way 

(*>) 

to define it, is to change in 5, (1.2.43), the exterior differential d to d , 
(2.3.14). The latter formula gives also the connection between the both 
operations, so that 



, (3) ( fc ) (3) _. (3) , (3). , 

div a := —* d * a = o a + * (r A *V T a ) (2.4.3) 



and finally 



j . (3) . (3) (3) 

div a:= 5 a -G» a . (2.4.4) 

The extension of this three-divergence operation to covectors which 
have not been beforehand projected onto three-space of the reference 
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frame, is achieved by a compensation of the redundant terms, as in 
the curl case. The final expression of the three-dimensional divergence 
operation in the monad language (the point of view of a reference frame) 
is 

div a := 5a — G • a — 2 a 6 — rd a . (2.4.5) 

It is worth mentioning that when the reference frame congruence rotates 
(the three-space is non-holonomic) , the second-order operation div curl 
does not vanish identically, but it yields 

(3) (3) 

div curl a = uj • £ T a . (2.4.6) 

This can be easily checked with help of relation (2.3.6) and the defi- 
nition of Lie derivative. We leave to the reader to compare this result 
with the identity (13.9) in Zel'manov and Agakov (1989); our next task 
is here to consider other important differential identities for accelera- 
tion, rotation, and rate-of-strain tensor, in particular those involving 
also curvature. 

The first of them is a scalar identity most easily obtained when one 
applies operator 5 to the rotation vector uj: 

Suj = uj a . a = 2uj»G. (2.4.7) 

This identity means that in the derivatives of uj which enter Su, the 
second derivatives of monad r cancel out. The same is true for the mag- 
netic field vector B which is a part of dA (A being the electromagnetic 
four-potential covector), and the second derivatives of A, cf. (4.3.11). 

Consider now curl applied to the acceleration of a reference frame. 
The resulting covector should lie in the three-space, hence all terms 
proportional to r must vanish. But first we have 

curl G = -2 * [G A *(r A uj)} - *(r A G) (2.4.8) 

where the first term is collinear to r, so that it will be exactly com- 
pensated, while the last term should be rearranged in such a way that 
it would not lead to a vicious circle. To this end we use here relations 
d(r AG) = -2dA (see (2.3.10)), 

£ T U = (uJ /3 . a T a + UJa^.^dx 13 



38 



CHAPTER 2. REFERENCE FRAMES' CALCULUS 



and 

Tq;/3 + T/3;a = G a Tp + G pT a + 2D af3 , 

from where 

icurl G = -{u a . a - 2G • uj)t + £ T u - 2u a D aP dx p + 2Quo. (2.4.9) 

Here the first right-hand side term (the parentheses) vanishes by virtue 
of (2.4.7), as we have expected. 

In order to find identities involving curvature, we already introduced 
alongside with the FW-curvature (1.2.33), also r-curvature making use 
of the operator (1.2.38) having remarkably better properties that those 
of (1.2.32). Thus 

R (w, v) := V;V; - V^V; - V\ wM . (2.4.10) 

This r-curvature has two pairs of skew indices, but there is no analogue 
of Ricci identities for it, so that one has to take a linear combination of 
its components in order to obtain a three-space tensor with properties 
of the standard Riemann-Christoffel tensor (cf. a different approach in 
Zel'manov and Agakov (1989)): 

1 M [r] [r] [r] [r] [r] 

r a p~,6 — g(2 Raf3~/5 +2 i? 7< 5a/3 + -R 7/ 3q<5 + Ra5~if3 + Ra~,f35 + R(36a-y)- 

Such a combination automatically satisfies Ricci type identities if the 

[r] 

building-material object (here, R) is skew in its first and last pairs of 
indices; this antisymmetry is also inherited by r a/3l s, yielding therefore 
symmetry under a collective permutation of [af3] and [yS], so that r^ 7 : = 
ra.p~f&b a6 (of the type of Ricci tensor) becomes symmetric. 

Since the operator (2.4.10) identically annuls r, we apply it to the 
basis vector b\ closely related to the monad. In contrast to r, we have 
here not a pure zero, but a non-trivial expression because of the "scalar" 
nature of the free index of b\ with respect to coordinate-free operations. 
Hence 
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R K)b x = [R a p jS b K AK b i + W + V)(A,a + A v \) 

-(D^ x + A, X )(D U K + A v K )]b K . (2.4.11) 

For the same reasons, the identity connected with r can be represented 

as 

W • [R(b x , b,)r] = R^r%l = 2®X - (£>/ + A^), a b«b{ - 2G p A px . 

(2.4.12) 

We call identities (2.4.11) the old generalized Gauss equations and 
identities (2.4.12) the generalized Codazzi equations (see Mitskievich, 
Yefremov and Nesterov (1985); for the standard theory of hypersur- 
faces which is not based on the reference frame approach, see Eisenhart 
(1926, 1972)); the generalization is meant in the sense that they are 
written in non-holonomic local three-submanifolds of rotating refer- 
ence frames. These equations turn into the usual Gauss and Codazzi 
equations correspondingly when u = (i.e., for a non-rotating con- 

M 

gruence we come to the usual theory of hypersurfaces). Since R does 
not possess the standard properties of curvature tensor, it is however 
necessary to make the next step of generalization using r a ^. Then the 
new generalized Gauss equations for the three-curvature read 

r K x^u = RawK^Wi + D Kfl D Xu - D KU D Xfl + A kX A^, (2.4.13) 

so that for the corresponding Ricci (three-)tensor we have 

r x , = r nXiXV b™ = R aPl 8b{blb aS + D;D Xu - D v u D XlI + A\A^ 

= R aPl sb{blb a& + D;D Xu - D v v D Xll - u x lu^ - u • u b x>1 , (2.4.14) 

where we have used simple identities A Kii A Xv — A RV A Xpu = A^ u A kX 
and A x A^y = cj a c^ + u • u b Xfl . 

The identities we have just obtained play a fundamental role in the 
3+1-splitting of Einsteins equations and in formulation of the Cauchy 
problem in general relativity; they are also important for a comparison 
of the latters with Maxwell's equations which we discuss in section 5.5. 
It is worth mentioning here another non-trivial fact: The definitions and 
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identities given above imply that in the case of a normal congruence 
(no rotation) there exists an integrating factor, i.e. such a function iV 
that d(N^ 1 r) =0. In a region with simple topological properties, the 
1-form ( = iV _1 r is not only closed but even an exact one, so that 
C = dt. If a quantity £ = Nt is simultaneously introduced (which 
usually appears as a contravariant vector), a simple reciprocal norm 
is realized, C ' £ = 1- Moreover, alongside with the trivial identity 
£^ = 0, a new important identity arises, £%Cn — 0- The projecting 
tensor (three- metric) possesses now the components bp = 5^ — £ Q C/3 and 
it is therefore constant relative to the Lie differentiation with respect 
to the vector field £: £^bg = 0. 

It is worth mentioning that similar ideas (although in spaces of other 
variables than the space-time ones) can be successfully formulated using 
Cartan's exterior forms, also in thermodynamics (Torres del Castillo 
1992) 



Chapter 3 



Equations of motion of test 
particles 



3.1 The electric field strength and mag- 
netic displacement vectors 

When the geodesic equation (1.3.29) is written in form 

u». a u a = 0, (3.1.1) 

its left-hand side is interpreted as an absolute derivative of the four- 
velocity with respect to the proper time of the particle (partial deriva- 
tive in special relativity if Cartesian coordinates are used). For a 
charged test particle the equation of motion can be written then in 
a general covariant form as 

mw%w Q = eF^Ua (3.1.2) 

(c/., e.g., Landau and Lifshitz (1971)). A comparison of this equation 
with the forms of geodesic equation (3.1.1) and (1.3.30) equivalent to 
each other, shows that the 3-form equation, 

u A *d(mu + eA) = 0, (3.1.3a) 

or its dual conjugate (1-form equation), 

*(u A *d(mu + eA)) = 0, (3.1.36) 
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are equivalent to (3.1.2). Here 



A : = 



a a dx^ 



(3.1.4) 



is 1-form of the electromagnetic four-potential, which should not be 
confused with the reference frame rotation tensor (2-form) A (2.3.5); 
anyhow, we shall use more frequently the rotation vector uj (see also 
(2.3.5)). 2-form of the electromagnetic field strength is then 



while m and e in (3.1.3) are the rest mass and electric charge of the 
test particle correspondingly, u being its four-velocity. 

Thus we come to a conclusion that the complete expression for the 
Lorentz force in the right-hand side of (3.1.2), can be represented as 



Consider now a question: how could one approach logically to the 
introduction of the concept of strength of a physical field? As one 
of the most simple cases studied in the theory of relativity in great 
detail, electrodynamics represents an ideal example of a formulation of 
such an approach. On its basis, one can insist that for a field acting 
on a corresponding test charge in some state of motion with a certain 
force, the field strength is a quantity which enters the expression of the 
force {pro unit charge value) and which does not include characteristics 
of the motion of the test charge. Thus the field strength describes a 
field existing in space-time independently of the charges used for its 
detection and measurement. In electrodynamics, this is the very tensor 
F; in a more complicated case of gravitational field, one has to introduce 
a concept of the relative field strength, since in the geodesic equation 
there is no gravitational force term whatsoever, so that it is necessary to 
consider the geodesic deviation equation (see (5.1.1)). Then the tensor 
rank of the relative strength becomes 4 and not 2, as it was the case 
for the field strength in electrodynamics. 

We know however from the special relativistic electrodynamics the 
connection existing between the four- and three-dimensional Lorentz 




(3.1.5) 



T = T^dx* = eF^ a u a dx^ = e * (u A *dA). 



(3.1.6) 
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force. This connection can be expressed from the point of view of the 
reference frame r as 

(3) It) 

T=u e{E + v x B), (3.1.7a) 

while 

(t) (t) 

F=ueE*v (3.1.76) 

(cf. the decomposition (2.2.5) and the three-projection of a vector, 
(2.2.4), as well as definition of the three- velocity of a particle, (2.2.11)). 
We rewrite now the first of these relations using the definition (2.2.15) 
of the vector product, 

t? E + *(u At AB) = *(t A *F) + *(v A *F) - [r ■ *(v A *F)]t, (3.1.8) 

where the right-hand side is also projected onto the three-space orthog- 
onal to r-congruence. This expression can also be written as 

E + v x B = F aP (T P + v 13 )^ . (3.1.9) 

This relation holds for particles moving arbitrarily with respect to the 
reference frame, i.e. for any vector v with a norm less than unity. Then 
one can write immediately (Mitskievich and Kalev 1975) 

E = F aP T P dx a (3.1.10) 

or equivalently 

E = *(t A *F) (3.1.11) 
(cf. (2.3.6)). From the arbitrariness of v the remaining part yields 

E^r^b" = F^bfr, (3.1.12) 

or better 

E KXlxv r K B x lf Ab u = -Fyjr A b v . (3.1.13) 

Here the left-hand side is equal to 2*{jAB). But since *(r A*(r AB)) = 
—B, we have after exterior multiplication of (3.1.13) by r from the left 
and a subsequent dual conjugation, 

B = *(tAF), (3.1.14) 



44 



CHAPTER 3. EQUATIONS OF MOTION 



an analogue of the expression (3.1.11) (c/. (2.3.5)). Then as an analogue 
of (3.1.10), we have expression with the opposite sign (Mitskievich and 
Kalev 1975): 

B = -F*^dx a . (3.1.15) 

Thus we have obtained very simple and applicable for all reference 
frames and arbitrary gravitational fields definitions of the electric field 
strength and magnetic displacement vectors, E and B. We remind 
that they follow unambiguously from the standard form for the Lorentz 
force, i.e. from the electric and magnetic fields interpretation in the 
spirit of traditional theory. The definitions (3.1.11) and (3.1.15) yield 
automatically a decomposition 

dA = F = E At + *(B At), (3.1.16) 

and vice versa (c/. (2.3.11)). 

In order to illustrate the simplicity of application of these defini- 
tions to concrete problems, we calculate now the electromagnetic field 
invariants through the three-space vectors E and B using the contracted 
crafty identities (1.2.6) and (1.2.7). To this end let us multiply each of 
the identities by t\t^ and make use of (3.1.10) and (3.1.15): 

F crT F crT = 1(B * B — E * _£7), (3.1.17) 

F* T F aT = AE • B, (3.1.18) 

where the definition of three-scalar product • (2.2.6), (2.2.8) is also 
taken into account. 

Let us build, with help of (3.1.16), self-dual quantities while first 
writing the expression dual conjugate to (3.1.16), 

*F = —B At + *(E A r). (3.1.19) 

This relation as well as its counterpart (3.1.16), can be more easily 
related to the component representation of the electric and magnetic 
fields (3.1.10) and (3.1.15) if one takes account of identities (2.2.20) and 
(2.2.22). 

We denote self-dual and anti-self-dual quantities using signs + or 
— , correspondingly, over the root letter: 

p:= F^fi* F =E Ar^fi* (% At) . (3.1.20) 
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It is easy to see that 



± 



(3.1.21) 



E:= E±iB. 



± 



The quantities F possess the properties 



± ± 
* F= ±i F ■ 



(3.1.22) 



Then the contracted crafty identities yield 



KvF^= -4i E • I, F^F^= 0, 



(3.1.23) 



the both relations following from (1.2.6) as well as from (1.2.7) since 
these mutually coincide in this case due to the equality (3.1.22). We 
see that an (anti-)self-dual electromagnetic field is equivalent to a si- 
multaneous existence of the both electric and magnetic fields, though 
differing one from the other by a complex factor (=p) , so that if for the 

initial 2-form F the correspondence F (E, B) holds, for F we have 



3.2 The energy and momentum of parti- 



We consider in this section general properties of the energy-momentum 
four-vector of mass points from the viewpoint of the theory of reference 
frames, as well as general problems of description of media consisting 
of such points which mutually interact via pressure and gravitational 



The general covariant energy-momentum vector of a mass point is 
decomposed into the energy scalar and momentum three-vector in a 
standard way (Mitskievich, Yefremov and Nesterov 1985) 




(3.1.24) 



cles 



field. 



p = mu = St + V, V ■ r = 0, 



(3.2.1) 



while 



S := p ■ t =u m, V :—P— £v, 



(3.2.2) 
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so that 

p-p = £ 2 -V»V = m 2 , £>m>0. (3.2.3) 

The scalar mass of a particle, often called its rest mass, is constant for 
elementary particles only (till a decay; moreover, if the lifetime of a 
particle is short, the constancy of m or universality of this quantity for 
different particles is limited by the uncertainty relation for time and 
energy. In the case of a macroscopic particle or a body, one cannot 
demand that its rest mass should be constant (e.g., for a radiating star 
or a body which throws away or evaporates its parts). Note that in 
general the idea of a fixed rest mass is alien to classical (non-quantum) 
physics, and it comes from observation of stable elementary particles 
for which this constancy has an exclusively quantum nature. 

If a medium consisting of particles which form a perfect (possibly, 
electrically charged) fluid, is considered, this can be characterized by 
the energy-momentum (or: stress-energy) tensor, 

T pf = (p + p)u® u-pg, (3.2.4) 

and the four-current density, 

j = pu (3.2.5) 

(not a vector density, but a vector!). Here p and p are invariant (taken 
in the co-moving reference frame of the medium) densities of the energy 
(mass) and charge of the fluid, and p is its pressure (not to be confused 
with the four-momentum in (3.2.1): we did not change the somewhat 
contradictory standard notations which in fact do not meet together 
in a close context). In such a co-moving frame the energy- momentum 
tensor takes the form 

T pf = pu <g> u - pw, (3.2.6) 

which justifies the choice of structure in (3.2.4); here the four- velocity 
u of the fluid plays the role of monad, and w = g — w ® -u is the 
corresponding three-metric tensor (signature of w is (0,-1,-1,-1)). 
Then in an arbitrary reference frame the energy-momentum tensor of 
a perfect fluid is 

T pf = { u p+Uu -Up r®r+ 
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+ u (/i + p)(r £g> v + v <g> r + v <g> v) — p b 



(3.2.7) 



where scalar coefficient in the first term represents the energy (mass) 
density in the reference frame r, the coefficient in the second term 
expresses the proportionality of density of the energy flow (the three- 
momentum density) of the fluid, to the three-velocity of the fluid in 
this reference frame; the last term describes the isotropic part of pres- 
sure (and stresses), while the next to the last one, characterizes its 
anisotropic part which is due to motion of the fluid with respect to the 
reference frame. It is interesting that invariant quantities, the energy 
density and the pressure, taken in the co-moving frame of the fluid, 
both contribute to the energy density, its flow density, and anisotropic 
part of the pressure, taken in an arbitrary reference frame, while the 
fluid undoubtedly is a Pascal one. By evaluating the involved quanti- 

(t) 

ties one has to take into account the fact that the relativistic factor u 
is always not smaller than unity (^u= dt/ds = (1 — v 2 )~ 1 / 2 > 1^. Ex- 
pression for the non-invariant densities of the charge and three-current 
are much simpler, 



It is worth emphasizing that the expressions "invariant" and "nonin- 
variant" are both used for scalar quantities depending not on the choice 
of a system of coordinates, but on the choice of a reference frame only, 
while speaking on "three-current" and "three- velocity", etc, we mean 
in fact four- vectors, which are lying in the local submanifold orthogonal 
to the monad r. One has to remember that no transformations of coor- 
dinates can influence the choice of reference frame, while in one and the 
same system of coordinates, there could be considered simultaneously 
as many different reference frames as one would desire. Exactly in this 
sense the energy density (even if non-invariant under transitions be- 
tween different choices of r) should always be a scalar under arbitrary 
four-dimensional coordinates transformations, though it might describe 
energy density of the fluid particles not in the state of rest. We speak 
here in such a detail on these elementary facts because they are in an 
acute contradiction with the vulgar and wide-spread understanding of 
reference frames as exclusively systems of coordinates. 



j —u p{r + v). 



(3.2.8) 
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Now let us touch upon the problem of dynamics of a fluid, which is 
expressed by vanishing of the four- divergence of its energy-momentum 
tensor, T pf/ ^ '.„ = 0. This means, of course, that we consider now an elec- 
trically neutral fluid, otherwise it would exchange energy-momentum 
with the electromagnetic field, so that we had to take the latter into 
account. Moreover, we shall consider here the fluid's behaviour only in 
a co-moving frame, u — t. From (3.2.6) we have 

/Vr Q = £ rf ji = -2(fji+p)e, (3.2.10) 
p, a b% = (gradp) = (fi + p)G . (3.2.11) 

We see that a non-homogeneity of the pressure produces a force acting 
on particles of the fluid and thus making their motion non-geodesic. 
But if the pressure is homogeneous in the co-moving frame of the fluid, 
its particles do move geodesically (c/. (Synge I960)), as it is the case for 
homogeneous cosmological models, otherwise the sum of the pressure 
and energy density should be equal to zero. We shall not give here any 
further discussion of these problems, and remark only that the last of 
the mentioned possibilities is closely connected with the properties of 
the de Sitter universe. 

3.3 Monad description of the motion of 
a test charged mass in gravitational 
and electromagnetic fields 

We consider here from the viewpoint of a reference frame r the split- 
ting of eq. (3.1.4) which describes the motion of a test charged mass 
point. Its non-gravitational part (3.1.5) has already met a comprehen- 
sive treatment in section 3.1 from where we shall borrow the corre- 
sponding expressions; now it remains to consider splitting of the term 

*(u A *dp) = u *[(r + v) A *d(Sr + V)] (3.3.1) 

in which it was suitable to pass to the dynamical quantities, the en- 
ergy 8 and three-momentum V characterizing the particle and already 
introduced in section 3.2. 
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It is easy to see that the four terms into which splits the right-hand 
side of (3.3.1), read 

*(r A *d(£r)) = d£- (£ t S)t - SG, (3.3.2) 

*(v A *d(Sr)) = 2£vxu-{£vG-v d£)r, (3.3.3) 
*(r A *dV) = -£ T V (3.3.4) 
(this is true only for a three- vector orthogonal to r, here, V), and 

*( v A *dP) = -(u • £ r P)r + - Pa, m )& A - (3.3.5) 

Differentiation of the identity (3.2.3) for m = constant, yields 

£ T (£ 2 ) = -£ T {g^V,V v ) = -1V»£ T Vp - V,V u £g^, (3.3.6) 

where £ T g^ = -(2D^ + r»G u + t u G") and v fM £ T V^ l = -v • £ T V, 
hence 

£ T £ = v £ T V + ED^vflv. (3.3.7) 

Note also that 

£ T £ — v • d£ — {m/£)£ u £. (3.3.8) 

As a result we obtain the splitted equations of motion of a particle 
having m = constant: two equivalent forms for the scalar equation (the 
component along r), 

£ T £ - v • d£ = £v^v u D^ + v(eE- SG) (3.3.9) 

and 

£ u £ = (£/m)[£ Vll v u D^ + v(eE- SG)], (3.3.10) 

and the vector equation (projection onto the three-space of the reference 
frame) , 

£ T V + V T V V = e(E + v x B) + £(-G + 2vx u), (3.3.11) 

cf. (Mitskievich, Yefremov and Nesterov 1985). In the last case, we 
do not employ the Lie derivative with respect to the four-velocity of 
the particle, since being applied to (co) vectors (in contrast to scalars), 
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this derivative brings the quantity out of the three-dimensional sub- 
space. In equation (3.3.11) the differentiation operation V£, (1.2.34), 
is used, which takes the reference frame explicitly into account. The 
corresponding term can also be written as V T V V = V KU v u b fi . 

In the resulting equations, a far-reaching analogy between gravita- 
tion and electromagnetism can be clearly traced in the aspect of the 
action of these fields on test particles. Below we shall further discuss 
this analogy (see chapters 4 and 5). 

The equations of motion (3.3.9) to (3.3.11) of course hold also in 
special relativity when no gravitational field (curvature) is present; then 
the terms containing G and oo represent the non-inertial frame effects. 
The very last term in (3.3.11) is then clearly the Coriolis force; the 
centrifugal force is hidden in the next to the last term proportional to 
G, if acceleration of the reference frame is expressed through oo and 
radius vector in the case of a pure (rigid) rotation. When gravitation 
is present, a mixture of gravitational and inertial forces arises, these 
forces being unseparable one from the other (remember the equivalence 
of gravity and acceleration, usually considered at a more naive level). 

3.4 Motion of photons, the red shift and 
Doppler effects 

Some concrete physical effects can be however calculated and described 
in such a way that they are expressed through scalar quantities being 
functions of states of motion of the participating particles only, in- 
cluding, if necessary, quanta which move with the fundamental speed. 
While the massive objects mentioned above may be used (if their world 
lines do not intersect) to build a reference frame in which these partici- 
pating objects are at rest, the massless quanta cannot of course be used 
in constructing such reference frames, since null lines are alien to the 
latters. Therefore one has not to unnecessarily universalize the trend 
towards description of all physical effects using reference frames. When 
they enter naturally into an experiment scheme or simplify essentially 
its description or calculations, they are surely advisable, but a mu- 
tual adjustment of setting an experiment and constructing a reference 
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frame, would be of an appallingly low standard, if this would be done 
for the end of realization of the idea of universal application of refer- 
ence frames only. As an example of such a situation we consider here 
a rigorous description of one of the most important relativistic effects, 
that of frequency shift of a signal, which in fact unifies the gravitational 
red shift and Doppler shift effects (in general, their unambiguous sep- 
aration is impossible). Some examples of alternative separation cases 
can be easily invented, but, if the space-time under consideration is 
non-stationary, there exists no criterion for choosing a preferential in- 
terpretation. Therefore we shall not be captivated by such a scholastic 
attitude, and we shall consider here an example of use of a local refer- 
ence frame, up to a numerical evaluation of the frequency shift effect. 
Such local reference frames should be chosen individually for every spe- 
cific problem, and there exists no general formalism of their application, 
in contrast with global reference frames for which the monad formalism 
excellently works. 

Four decades ago, Schrodinger (1956) has proposed an invariant 
procedure for determination of the frequency shift effect (in the de- 
scription of the cosmological red shift), but this procedure remained 
practically unnoticed until Brill (1972) has given in his short report a 
rigorous deduction of Schrodinger's formula using coordinates-free ap- 
proach, the same as applied in our paper. While giving here a deduction 
of the Schrodinger-Brill formula, we propose to the reader to draw by 
her/himself the corresponding figure, so simple the idea is and so ad- 
vantageous this would be for its best understanding. Let the emitter 
and detector of signals have time-like world lines 1 and 2 correspond- 
ingly, while the signals propagate along null geodesies from the emitter 
to detector (thus the light cone has its apex on emitter's world line, and 
the generatrix intersects the detectors's world line). Then two signals 
sent from the emitter with an interval of the proper time between them, 
equal to some standard period, reach the detector being divided by an- 
other interval of the proper time (now, along detector's world line), 
different from the initial standard period. If the detector is provided 
with a frequency standard moving together with it and reproducing 
the standard period, now from the viewpoint of the proper time of the 
detector, it is possible to compare these two time intervals, i. e to mea- 
sure the frequency shift by the transmission of the signal from emitter 
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to detector. This is a fundamental principle of the (general as well as 
special) relativity theory that the period (or, of course, frequency) of 
signals expressed through the proper time, is one and the same, inde- 
pendent of the choice of the world line of the generator of these signals 
on which the proper time is to be measured. By the way, the same 
idea of relativity (and the equivalence principle), makes it in principle 
impossible to separate Doppler shift and gravitational red sift effects. 
By definition, the red shift is expressed via the formula 

Z = (l>em - ^dcOMm = 1 ~ ds 1 /ds 2 (3.4.1) 

where for a computational convenience, we introduced infinitesimal 
proper time intervals on the corresponding world lines, 1 for the emit- 
ter and 2 for detector, while ds 2 is the interval between two signals 
which reach the detector, having initially — at their departure mo- 
ments on the world line of the emitter — some standard proper time 
interval ds\ between them. It is clear that this formula simultaneously 
describes the both usual gravitational red shift effect and the Doppler 
shift (due to motion of the emitter as well as the detector). In general 
one cannot speak on a relative velocity between detector and emitter, 
if the space-time is curved: there will be no unambiguous way to com- 
pare the four- velocities since the comparison would depend crucially on 
the choice of the transportation path, so that only in the case of the 
Minkowski space-time (absence of any genuine gravitational field) such 
an absolute comparison of motions is possible: this is why we are so 
much accustomed to the independent existence of the Doppler shift in 
special relativity. Thus our task now is to determine the connection 
between ds\ and ds 2 for one and the same pair of signals. 

Let us complement the time-like world lines 1 and 2 by introducing 
some intermediary time-like lines (the concrete choice will not affect 
the results) which, as the previous lines 1 and 2, have not necessarily 
to be geodesies. Thus we have built a family of lines which may be 
parametrized holonomically when a parameter a will vary along each 
of them (the other parameter, A, is introduced a little later, and it 
enumerates the individual lines continuously). Then the tangent vector 
field of this family of lines is q = d a . The other family of lines will be 
that of the null lines of the signals propagating between the initial world 
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lines 1 and 2. These null lines are, of course, geodesies, with a canonical 
parameter A and tangent vector field k = d\, so that 

V k k = 0, k-k = 0. (3.4.2) 

Let the parameter A be holonomic too, thus the condition to have 
two holonomic parameters (which may be used as coordinates) can be 
expressed as 

V q k - V k q = [q, k] = (3.4.3) 

(cf. the sixth axiom of the covariant differentiation). Note, on the one 
hand, that the scalar product q ■ k is constant under transport along 
the null lines. In order to prove this fact, one has to differentiate q ■ k 
with respect to the parameter A: 

d x (q ■ k) = V fc (g • k) = k ■ V k q = k ■ V q k = ^7 q (k ■ k) = 0. (3.4.4) 

On the other hand, this is not the vector field q (introduced for the first 
family of lines) which is physically meaningful, but the four-velocity 
vectors u on the world lines 1 and 2, or equivalently, the proper time 
intervals along these lines. Therefore it is worth writing the conserved 
(by the null transport between the lines 1 and 2) quantity as 

q ■ k = (ds/da)u ■ k. (3.4.5) 

The holonomicity property of a, A means in particular that the dif- 
ference da is one and the same on the both lines 1 and 2, if it is 
taken for one and the same pair of signals; then evidently (u ■ k ds)i = 
(u ■ k ds) 2 . We have thus finally obtained the Schrodinger-Brill formula 
(Schrodinger 1956; Brill 1972; Mitskievich and Nesterov 1991) 

z=l-(u-k) 2 /(u-k) 1 (3.4.6) 

which holds both in special and general relativity and describes Doppler 
frequency shift in the Minkowski world, as well as the gravitational red 
shift effect in the Schwarzschild space-time, and the cosmological red 
shift in the Friedmann world, to mention a few of the principal cases. 
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We now apply this general formula to two cases, the cosmological 
red shift and a free fall of the detector onto a black hole. In the first 
case the metric is 



It is most simple to use the coordinates co-moving with the matter 
o universe (with which we also connect the both emitter and detec- 
tor of radiation). Let the detector rest in the spatial origin (x2 = 0; 
this is justified by homogeneity of the Friedmann universe), so that 
signals propagate radially. Their equation of motion, ds 2 = 0, is inte- 
grated trivially: r\ 2 — r\\ + Xi- Here rji and r/ 2 are cosmological times 
of transmission and reception of a signal, the difference between them 
being equal to the cosmological coordinate "distance" of the radiation 
source from the origin where the observer resides (this xi is & t the 
same time integration constant in full accordance with the use of a co- 
moving frame in which all the matter of the universe is at rest). Then 
u = a~ l d v , k = di] + dx (it is easy to check that k forms a geodesic 
field). A substitution into (3.4.6) yields z = a 2 j 'a\. If expansion of 
a (771) about the observation moment r\ 2 is performed (with the corre- 
sponding assumption that the two moments are sufficiently close one 
to the other), we obtain a\ = a 2 — (da/dr)) 2 Xi- Since the physical dis- 
tance between the observer and radiation source is equal to / = ax, this 
formula yields the standard Hubble expression 



(where transition to the physical time of the observer is also performed 
according to the relation adr] = dt). In this example a system of coor- 
dinates (co-moving with the reference frame whose language is used in 
our discussion) was employed in which the both radiation source and 
observer were continuously at rest, moving at the same time geodesi- 
cally. Nevertheless, the obtained effect is treated usually as a Doppler 
shift. It is closely connected with the velocity of galaxies' parting away 
(the expansion of the universe); however, in our frame the matter is 
everywhere at rest, but the scales are changing (the rate-of-strain ten- 
sor is non-zero for the reference frame under consideration). But it is 
clear that there can be no physical sense in arguing into one or another 




(3.4.7) 



H = a 2 da/dr] = a 1 da/dt 



(3.4.8) 



3.4. RED SHIFT AND DOPPLER EFFECTS 



55 



interpretation of the nature of this red shift effect, the gravitational or 
Doppler one. 

The second case is related to a Gedankenexperiment concerning a 
situation when over a black hole, a mother craft is hovering, and a probe 
dives from it with zero initial velocity into a free fall, while standard 
signals are continuously sent to it from the mother craft (Mitskievich 
and Nesterov 1991). The problem is to find out the frequency shift of 
these signals received by the probe. We consider the general case of a 
black hole when it possesses a mass M, Kerr parameter a, and electric 
charge Q. Determine first of all the tangent vectors u and k in the 
Boyer-Lindquist coordinates. The first integrals of geodesic equation 
in the Kerr-Newman field are (see Mitskievich and Nesterov (1991)) 

ASi = EA, 

£ 2 r 2 = E 2 {r 2 + a 2 ) - A(K + W 2 ) =: F 2 , 

T?9 2 = K -r]a 2 cos 2 6 - E 2 a 2 sin 2 6, 

AS0 = aE(2Mr - Q 2 ) 

where A = r 2 + a 2 + Q 2 - 2Mr, A = (r 2 + a 2 ) 2 - Aa 2 sin 2 0, £ = 
r 2 + a 2 cos 2 8, E and K being the first integrals of energy and the 
total angular momentum (determined with help of the Killing tensor, 
cf. (Marck 1983)), i] = or 1 depending on null or time-like lines are 
considered; the speed of light and Newtonian gravitational constant 
are assumed to be equal to unity. Let now 9 = 0, so that the motion 
is radial up to such dragging effects (in <fi direction) which cannot be 
excluded at all radii simultaneously. 

Although the world line 1 is non-geodesic, we choose it so that it 
corresponds to all other motions, i.e. u\ = i d t + 4> d^. Here the expres- 
sions from above are taken for t and <fi while i] = 1, while normalization 
of the four-velocity yields a special choice of the constant E (the Carter 
integral K is already fixed by the demand = 0). The same value of E 
is inherited by the probe since the world line 1 gives in fact the initial 
data of its motion. For the probe we have then u 2 = t d t + r d r + <j) d^, 
with the same remarks which were made with respect to u\. Since 
the angle 9 does not change along the world lines, the coordinate r 
is essential only (all quantities do not explicitly depend on t and 0), 



56 



CHAPTER 3. EQUATIONS OF MOTION 



so that it is not necessary to determine the intersection points of the 
lines under consideration. In the common (initial) point of the world 
lines 1 and 2 (only the latter is a geodesic) the both lines have a com- 
mon tangent vector. As to the null geodesies, their tangent vectors are 
k = t d t + r d r + (pod^, and since all components are proportional to 
one and the same constant E , this can be equalized without loss of 
generality to the already used E. Hence to = t, 0o = <fi, 9 = 9 = 0, 
and the only remaining non-trivial relation is 

£ 2 f 2 = E 2 [(r 2 + a 2 ) 2 -Aa 2 sin 2 6}. 

Thus SV 2 , = S 2 r 2 + AS, and k = u<i + (ro — f)d r . The radius at which 
the signal comes to the world line 2, we denote as the moving coordinate 
r. Since r\ = 0, [k ■ u)\ = 1, and in general k • u — 1 + r(r — r)g rr . It 
is easy to find that r = — F/S and r = — EA 1 / 2 /Y,, where the choice 
of sign (minus) is determined by the direction of motion (the fall). 
Moreover, since AE 2 = F 2 + AS, we obtain for the red shift measured 
by the probe falling onto the Kerr-Newman black hole, finally 

z = (F 2 /SA)[(1 - SA/F 2 ) 1/2 ]. (3.4.9) 

Since at the horizons (the both event and Cauchy ones) A = 0, we 
see that the limiting value of the red shift is z — | when the probe 
crosses a horizon. It is remarkable that this quantity does not depend 
on anything whatsoever — on mass, Kerr parameter and charge of the 
black hole, and on the radius at which the mother craft is hovering (with 
a drift in the azimuthal direction) and the polar angle value to which all 
the motions correspond. The frequency shift occurs namely to the red, 
i.e. it can be treated as if the Doppler shift would dominate over the 
gravitational (in this case, violet) one. Thus, using the red frequency 
shift of the standard signals coming from the mother craft, it is possible 
on the probe to determine the moment of crossing the horizon of any 
black hole. In this example, in contrast with the cosmological case, it 
is impossible to introduce a frame which is simultaneously co-moving 
with respect to the both source and receiver of the signals: at the initial 
moment their world lines coincide, so they cannot belong to one and the 
same congruence. This is however a meaningful example of application 
of a bilocal approach to reference frames, complementary with respect 
to the monad formalism. 
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3.5 The dragging phenomenon 

The dragging phenomenon embraces a large variety of effects predicted 
for test particles moving in gravitational fields of rotating (and simi- 
lar to them) systems of material objects, and it probably has a very 
profound and general physical nature. It is usually characterized as 
"dragging of local inertial frames" , being considered — in the absolute 
majority of cases — in space-times admitting a time-like Killing vector 
field £ (see, e.g., Misner, Thorne and Wheeler (1973)). It is clear that 
in such space-times a privileged reference frame is easily introduced, 
namely that which coincides with the Killing congruence. Rotation of 
this reference frame describes then the dragging effects. We shall see 
that this rotation (essential stationarity of the space-time) is a neces- 
sary prerequisite for existence of dragging, though usually in discussions 
of the dragging phenomenon nothing is mentioned concerning reference 
frames being used, maybe some traditional words about the local refer- 
ence frames only (but without any elements of a sensible formalization 
of the reference frames theory). In this section, we attempt to give an 
introduction to consideration of the dragging phenomenon in the realm 
of this theory, as well as additional hints pointing at some prospects 
and needs in generalizing the very concept of dragging. 

In a space-time admitting time-like Killing vector, there exists a 
preferred reference frame with 

" (3.5.1) 



Such a Killingian reference frame can be interpreted as a co-moving one 
with the gravitational field (or with the physical system itself which pro- 
duces this field), so that the frame is in the state of acceleration and 
rotation together with the space-time being considered (no deforma- 
tion can obviously be present in this case). If this space-time is not 
static, but stationary one, the dragging phenomenon occurs. It is con- 
nected with existence in such a reference frame, of a gravitomagnetic 
field — the rotation vector uo (cf. Nordtvedt (1988); Jantzen, Carini 
and Bini 1992). If we change to another monad, f, this vector field 
may be of course transformed away, but since the preferred Killingian 
frame is supposed to be in a rotation, uj has to be considered as pos- 
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sessing an absolute meaning. We speak then on a rotating Killingian 
reference frame, whose rotation yields gravitomagnetic effects, and this 
is precisely the dragging phenomenon. Since there usually exists an- 
other (now, space-like) Killing vector field rj which may be considered 
as complementary to £, a combination of these two fields (with constant 
coefficients) also represents a Killing vector field, the concrete choice 
of coefficients determining a hypersurface on which the corresponding 
rotation vector vanishes locally (this analysis may be most character- 
istic for studies of the ergosphere region, e.g., of the Kerr space-time, 
and it corresponds to the local stationarity property of that space-time; 
cf. Hawking and Ellis (1973), p. 167. In the case of a pencil-of-light 
space-time, local stationarity is in fact a global property). 

It is worth giving here a concrete, but general enough example of 
how works the existence (or absence) of the rotation vector, from the 
point of view of a coordinated basis connected with r. The Killing 
equation (1.2.42) can be rewritten as 



One may choose coordinate (say, t) lines coinciding with the inte- 
gral curves of the field £, £ • £ > 0, together with a corresponding 
parametrization along these lines, so that in the new system of coordi- 
nates, £ = dt = St 8^. Then the last equation takes form g^t = (i.e., 
all metric coefficients are independent of the Killingian coordinate t). If 
there exist several independent Killing vector fields, the independence 
of all components g^ u of all corresponding new (Killingian) coordinates 
may be achieved, only if these Killing vectors mutually commute. If 
not, the constancy of metric coefficients is realized in different systems 
of coordinates only (not unifiable into one and the same system). We 
confine ourselves to a consideration of one single Killing vector field, £, 
and we will show that a non-rotating Killingian congruence corresponds 
to a static space-time, and vice versa. In the corresponding Killingian 
system of coordinates, 

ds 2 = gudt 2 + 2g ti dtdx l + g^dx % dx\ 

i,j — 1, 2, 3 if t — x°; all g^ are independent of t. Equivalently, 




(3.5.2) 




(3.5.3) 
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or 



ds 2 = {r^dx^) 2 + b^dx^dx' 



where 




9tn 



(The first of these last expressions represents a monad in coordinates 
co-moving with the reference frame (Zelmanov's chronometric invari- 
ants formalism), while the second one simply duplicates (2.2.1).) A 
complete separation of t and x % in (3.5.3) — orthogonalization of t axis 
with respect to all other (here, spatial) coordinates — occurs when the 
expression ( := r/^/g^ = dt + (gu/gtt)dx l is a total differential. This 
case is known as the static field case (a simultaneous orthogonality 
of t with respect to other axes and t-independence of all g^). This 
corresponds to fulfilment of the condition 



the same as vanishing of (or, equivalently, of a;). If u ^ 0, the 
space-time is a stationary (but not static) one. 

Thus we have come to the following chain of conclusions: It is always 
possible to choose the coordinates in such a way that £ = d t ; then 
£ • £ = 9tt, so that r = f / ^/gj t . Now, C := r/ yTgJt = {/ 9tt = Oj 
but it is possible to introduce a new time coordinate, t, so that ( = dt, 
if and only ifuj — 0. 

We shall consider first a very simple case of dragging, that which 
occurs in the equatorial plane of the Kerr space-time. The Kerr metric 
in the Boyer-Lindquist coordinates reads 




* 2 = (> " ? 



2 / ymr 




r 2 + a 2 cos 2 $ 



-dr 2 -(r 2 +a 2 cos 2 §)dd 2 



2 + a? cos 2 d 



r 2 — 2 r ymr + a 




r 2 + a 2 + 



2 r yma 2 r sin 2 •& 



r 2 + a 2 cos 2 $ 



) 



sin 2 tfd0 2 + 2 



r 2 + a 2 cos 2 1? 



2^mar sin 2 



d0rft, (3.5.4) 



so that 
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For a circular equatorial orbit, radial component of the geodesic equa- 
tion 

d ( dx»\ 1 dx a dxP 

d-s{ 9r ^) =0 = 2 9a ^^ (3 - 5 - 5) 

gives enough information for solving the problem. Here dx^/ds = 
(dt/ds)5t + (d(p/ds)S^ (r and d do not change in the course of mo- 
tion). There are two roots of the eq. (3.5.5), 




(3.5.6) 

V V ' T" L J 
(<j) := d(j)/dt), so that 

T± := 2tt I 4>± \~ x = 2ir (J— ± a] = T N ± AT (3.5.7) 



7m 



(we consider here a being a smaller term than ^/r 3 / (7m) in a physi- 
cally meaningful region of the motion). The main part of the period has 
here the standard Newtonian value T/v, while the dragging term sim- 
ply simply is additive. It is rather remarkable that it does not depend 
on the orbit's radius, mass of the attracting centre, and even on the 
gravitational constant, thus being suggestible from purely dimensional 
considerations. At the same time, it is worth stressing that this result 
is exact, and not approximate one (see Mitskievich and Pulido (1970), 
Mitskievich (1976), Mitskievich (1990); for an approximate approach, 
see Mitskievich (1979)). The fact that we have used here the coordi- 
nate (not proper) time and (coordinate) angle 0, does not make our 
results non-co variant, since these coordinates are the Killingian ones, 
thus invariantly reflecting geometric properties of the Kerr space-time. 
It is clear that the middle point between the two test particles moving 
in the opposite directions, then itself rotates with the angular velocity 

9m = 5 ~ (3.5.8) 

7mcr — r 6 



(in (Mitskievich and Pulido 1970), this point was interpreted as "the 
meeting point" of the particles), and a simple calculation shows that 



3.5. THE DRAGGING PHENOMENON 



61 



(Mitskievich 1976) if for usual celestial bodies the effect is quite small, 
for super-dense objects (such as pulsars) it is very large: see the Table 
in which we have taken the circular orbit's radius equal to r = "S^^R 
where R is radius of the real attracting centre, D the meeting point 
drift in seconds of arc per century (the Schwarzschild precession A in 
the same units being also given as illustrative data, the quantity r being 
then the large semi-axis of an elliptic orbit). 



Central 
body 


m, g 


R, cm 


L, 


g-cm l 

s 


T N , s 


D 


A 


Sun 


2 • 10 33 


7- 10 10 


1.3 


• 10 49 


2 • 10 4 


600 


10 6 


Earth 


6 • 10 27 


6 • 10 8 


7- 


10 40 


9- 10 3 


4.4 


670 


Jupiter 


2 • 10 30 


7- 10 9 


7- 


10 45 


2 • 10 4 


280 


10 4 


Rapidly 
















rotating 
















class B star 


3 • 10 34 


3.5 • 10 11 


2- 


10 53 


5- 10 4 


7.6 • 10 4 


10 6 


Pulsar at a 
















breaking 


10 33 


10 6 


4- 


10 48 


10" 3 


93 ^ 


750 sis' 




sec 


sec 


point 

















A closely analogous, but exotic dragging effect is present in the NUT 
field — the field of a gravitational dyon (Mitskievich 1981, Mitskievich 
1983). Here a gravitoelectric field of the Schwarzschild type is accom- 
panied by a gravitomagnetic field, similar to the magnetic field of a 
magnetic monopole. If one considers in this field circular orbits of test 
masses, one finds that the field centre cannot be in the plane of such an 
orbit, but it is shifted in the direction perpendicular to the plane, by 
a distance proportional to the value of the orbital angular momentum 
and to the NUT parameter /, the direction of the shift being the same 
as direction of the angular momentum L for I > 0, and opposite to this 
direction for I < 0. The effect is described by the exact formula 

costf = ~2Ey. 

Another analogous dragging effect is that of a Kerr gravitational 
lens when two light rays passing by the Kerr centre in its equatorial 
plane with the same impact parameter but on opposite sides of the 
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centre, are deviated differently by the gravitational field of the rotating 
source, and they are focused off the "straight" line going through the 
centre parallel to the initial common direction of the rays (a transverse 
shift of the focus of these rays). The magnitude of this shift from 
the above "straight" line is equal to the Kerr parameter (if measured 
in the units of distance), this being however an approximate result 
which becomes more precise for large values of the impact parameter 
(Mitskievich and Gupta 1980, Mitskievich and Uldin 1983, Mitskievich 
and Cindra 1988). 

Now let us consider a stationary case of the pencil-of-light space- 
time: 

ds 2 = dt 2 - dp 2 - dz 2 - p 2 d<p 2 + ^e\n(ap)(dt - dz) 2 . (3.5.9) 

Here e is the energy (or, equivalently, momentum ^-component) linear 
density, while a serves simply for making the argument of logarithm 
dimensionless (in a non-stationary case, both e and a are treated as 
arbitrary functions of the retarded time t — z). 

The already mentioned stationarity property of the metric (3.5.9) is 
non-trivial. There always exist Killing vectors £ = d/dt and rj = d/dz 
(as well as d/d(p, of course). The squares of £ and rj are not sign- 
defenite, so that t is not always a timelike coordinate, and z not always 
a spacelike one. In the region —1 < 87elno-p < +1 the coordinates 
are respectively time- and spacelike. When 876 map < — 1 the both 
coordinates are spacelike, and when 87eln<rp > +1, they are timelike. 
But at any point one may take linear combinations (with constant 
coefficients) of these two Killing vectors, yielding two new Killing vec- 
tors being correspondingly timelike and spacelike. Moreover, the band 
— 1 < 87elnerp < +1 may be then shifted to every position (retaining 
its finite width) if the following coordinate transformation is used: 

t = (1 - L)t' + Lz', z = —Lt' + (1 + L)z', L = -4ryeln(a'/(r), 

while 

87eln(<rp) — > 87eln(a / p)- 



From the geodesic equation it is easy to conclude about existence 
in the pencil-of-light field of the phenomenon of dragging; at the same 
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time, the timelike Killing vector congruence (in the above-mentioned 
band) rotates, so that in the privileged frame there exists a gravitomag- 
netic field (alongside with a gravitoelectric one). With A as a canonical 
parameter, the first integrals of motion are 

(1 + 8 7 eln(ap))^ = a, (1 - ^el^ap))^ = (3, 

2 d(j) 

p Tx=^ 

a, (3 and \i being the integration constants. For p we have the equation 

d 2 p 4 7 e 2 fi 2 

dX 2 p p 2 

(it can be easily integrated to a first integral, but this reduces to the 
squared interval (3.5.9)). In the right-hand side, the first term describes 
gravitoelectric attraction to z axis, and the second one, the centrifu- 
gal repulsion. The dragging effect is in this case compensated by the 
(conserved) ^-component of linear momentum of the test particle. 

If the test particle even did momentarily move without changing its 
^-coordinate, being however not on a circular orbit in a plane perpen- 
dicular to z axis, this particle is dragged in the next moment along this 
axis, the sign of the dragging being determined by the expression 

(a - /3)8 7 e^ 
Po 

where po is the value of p when, momentarily, dz/dX = 0, and Ap is 
the deviation of p from p . 

It is worth mentioning that a photon (lightlike geodesic) moving par- 
allel to z axis in its positive direction, does not interact with the pencil- 
of-light field at all (see for a general discussion Mitskievich (1989)), but 
a photon moving in the opposite sense will fall subsequently onto the 
pencil of light (z axis). This fact of no interaction holds not only 
for test objects in the pencil-of-light field, but it is also true for self- 
consistent problems with light-like moving objects in general relativity. 
Thus we come to exact additivity of gravitational fields of parallel (not 
anti-parallel!) pencils of light; this being an exact generalization of the 
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known approximate result obtained by Tolman, Ehrenfest and Podolsky 
(1931). 

Another case of the pencil-of-light field involves not only luminal 
(null) motion of the field's source, but its rotation too, this rotation 
being considered as an analogue of polarization of an infinitesimally 
thin beam of light. The corresponding metric reads 

ds 2 = 2d(t - z)[d(t + z) + (47eln(crp) - a<j>) d(t - z)\ - dp 2 - p 2 d(p 2 

(3.5.10) 

which may be modified to 

ds 2 = 2d(t - z) [d(t + z)+ 47e \n(ap)d(t -z) + (t- z)ad(f>] - dp 2 - p 2 d<f) 2 

(3.5.11) 

where a is a new constant corresponding to the "polarization" (in a 
more general case considered from another point of view in chapter 
5, instead of the product (t — z)a there stands an arbitrary function 
g(t — z); then the space-time becomes however non- stationary). The 
metric (3.5.10) describes a stationary gravitational field, thus permit- 
ting application of the standard approach to the dragging phenomenon. 
Here, as well as in the alternative version (3.5.11), we see that the ro- 
tational motion when coexisting with the luminal one, inevitably leads 
to a dragging of dragging effect, that is, to appearance of the prod- 
uct term dzd(f), in addition to dtd<f>, dtdz being already present by the 
virtue of the pencil-of-light nature of the initial metric. This is not a 
merely accidental presence of all three terms in a concrete solution, but 
inevitability of the presence of them all, if any two of them are already 
present in a solution. 

The simplest approach to consideration of dragging in the pencil-of- 
light space-time consists of a calculation of the (proper) time derivatives 
of coordinates of a free test object. We mean here not covariant, but 
simply partial derivatives which manifest the tendency of varying the 
coordinates (or their derivatives, if only the corresponding higher order 
derivatives become non- vanishing) . Thus, considering as an initial con- 
dition an instantaneous state of "rest" , we see that the geodesic equa- 
tion yields (through the second order derivatives) a (non-covariant) 
acceleration towards the pencil-of-light source (in the — p direction); 
these (and possibly the third order) derivatives show, moreover, that 
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a tendency of dragging exists, acting in the <fi (polarization) and z 
(orientation of the light-like motion of the pencil of light) directions. 
We encounter then, as manifestations of dragging in this special kind 
of field, the germs of acceleration (see Mitskievich and Kumaradtya 
(1989)). 

Let us discuss this problem quantitatively using a more general form 
of the spinning pencil of light (SPL) metric, 

ds 2 = 2dv (du + k(v) \napdv + g(v)d(f) — dp 2 — p 2 d<j) 2 

(a possible dependence a(v) is inessential); here the variables v and u 
correspond to t — z and t + z in (3.5.10)). 

Since the metric coefficients are independent of u anf <fi, one imme- 
diately arrives at two first integrals of motion, 

dv dv d(f> 

9vu Tx =a>0, gv +- +g „- = -p. 

The further exact integration of the geodesic equation is in general 
impossible, but we are not interested in obtaining approximate solutions 
(e.g., using perturbative procedures), so another way to deal with the 
dragging phenomenon should be chosen. We shall consider only germs 
- the tendency of test particles motion, i.e., the first non- vanishing 
higher order derivatives of the particles' spatial coordinates when the 
initial state of motion is given. This initial motion is most naturally 
chosen as a state of momentary rest (if the test particle has non-zero 
rest mass, thus d\ = ds). Then 






We may now express z in analogy with the Peres wave: z = (l/2)(u — v) 
(cf. also (3.5.10)). Then 





o \ / o 

The remaining components of the geodesic equation yield 

/ d 2 p\ ka 2 
\ds 2 ) p 



(3.5.12) 
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a2Unap (3.5.13) 



(3.5.14) 

o r 

where the right-hand side quantities are also taken at the initial moment 
(in the sense of (...)o)- Since k = 8jeo > (the linear energy density 
of the SPL source must be positive), (3.5.12) describes attraction of 
the test particle to the SPL. The noncovariant accelerations in the 
directions of z and 0, (5.3.13) and (5.3.14), may be either positive or 
negative depending on values and signs of g, g, and k, as well as on p. 
The sign of dragging in the direction coincides with that of g, and this 
means that not the function g, but its first derivative is directly related 
to the angular momentum of the SPL which is responsible for dragging 
in this direction. Moreover, if g is constant, it can be transformed away 
by merely introducing u = u + g<f>, so that only g can have a physical 
meaning. 

As to the motion of a lightlike particle, its initial state should be 
now chosen in a different way, since the former initial conditions are in- 
compatible with ds 2 = 0. So we shall take only two of these conditions, 





Inserting them into ds , we have 



o 

which leads to two possibilities: 



dv \ / du ^ dv 
d\ ) „ \d\ ^ dX 





^ =0 



and 

(s + »»"5) - a 

The first one corresponds to absence of interaction between parallelly 
moving lightlike objects (in fact, also not test ones). The alternative 
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case (3.5.15) does not admit (dv/d\)o = since such would 
correspond to a world point (an event) and not a world line. Combining 
the condition (3.5.15) with the geodesic equation (including the first 
integrals of motion), we come to the relations 

/ d 2 p\ ka 2 



d\ 2 J Q p 

„2 




o 

2 



P 2 



d 2 u\ gga 2 , , 

d\ 2 =-—- ka lnap - 



o 

These accelerations are in fact the same as in the case of a massive 
test particle, though we have here u instead of z. We see that a photon 
moving parallel to a SPL does not interact with it (in particular, it does 
not feel dragging in the <fi direction), while a photon moving antiparallel 
to it, both falls onto the SPL and starts to rotate in the direction. 

Let us now apply definition of the rotation covector uo (2.3.5) to the 
Killingian frame in the Kerr space-time. In fact, this is the only object 
which could be connected with the quantitative manifestations of the 
dragging phenomenon. Since in a stationary space-time (as also in a 
static one) £ = d t , (3.5.1) reads 



9,a -_dx*; 



thus (2.3.5) yields 



9nt 



g^xE^ v dx v . 



First, we observe that k, /i — t, and A = r, i?. Hence, v — r 
(complementary to A). Then it is clear that on the equator, u is directed 
along i.e. the rotation occurs in the (±) <fi direction. This is quite 
natural and corresponds to a similar rotation of the source of Kerr's 
field. The general exact expression is actually 



u = -E tr ^ 
2 



9t<j>,0 — ^-gtt,v \ 9 T — \ Cjt<j>,: 




(3.5.16) 
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(contravariant vector representation). 

On the equator (i? = n/2), the rotation vector (3.5.16) takes the 
form 



We have to compare this absolute value (which is a scalar) with (3.5.8) 
(which is neither a scalar nor component of a vector). They coincide 
(up to the sign) when the field is weak (i.e. r ^> 7m, r ^> a). Why this 
coincidence realizes for a weak field only? The answer could be that 
(1) we compare quantities with different geometric natures (see our 
remarks in the parentheses above) which may become similar in the 
vicinity of flat space-time only; (2) the very description of the dragging 
effect as a mean value of two rotation frequencies is highly tentative, 
and could serve well in the weak field approximation, but not in the 
general case. 

In this comparatively short account on dragging we can but mention 
many other aspects of dragging, though they can be easily related to 
the reference frames' theory. These include, first of all, effects of change 
of orientation of test bodies possessing multipole (beginning with the 
dipole) moments (thus without spherical symmetry, of scalar property). 
This is, in particular, the Schiff effect (precession of a gyroscope). Ef- 
fects of this type are sometimes classified as non-dragging ones, since 
they are connected with a local rotation of test particles without neces- 
sarily a translational motion, but we cannot consider this distinction as 
essential. See, e.g., Sakina and Chiba (1980), Marck (1983), Tsoubelis, 
Economou and Stoghianidis (1987), Tsoubelis and Economou (1988). 
This motion includes precession-like manifestations not only of drag- 
ging, but also of the tidal forces etc (see also Ashby and Shahid-Saless 
(1990)). The mentioned papers contain additional bibliography on the 
subject. 

A specific character of these effects includes arbitrariness of the gy- 
roscope world line which can be assigned artificially (usually, one takes 
a geodesic which is a good approximation for a freely moving gyro- 
scope). As another example one may consider a gyroscope suspended 
at some fixed point over the surface of the Earth. 



cq 




(3.5.17) 
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To this type of dragging belong effects involving non-scalar (in par- 
ticular, rotating) test bodies which do not change their orientation (rel- 
ative to the directions determined by symmetries of the gravitational 
field) in course of the motion. Thus no precession-type effects are now 
considered. For example, a spinning particle can hover over a Kerr 
source pole if 



. w 9 9x ^ „3z 2 — a 2 / 2m\z\ 

M(z 2 - a 2 ) = -2aS— 9 T 1 + 

z l + a 1 \ z l + a- 




Here z is position of the test particle over the pole and S, spin of 
the particle. This effect represents the spin-spin interaction (Epikhin, 
Pulido and Mitskievich 1972). An example of the spin-orbital interac- 
tion can be given for a circular motion in the Schwarzschild field of a 
pair of spinning test particles with antiparallel spins orthogonal to the 
plane of their common orbit (see the same publication). Their periods 
of revolution are approximately given by the formula 



/ / r 3 35 1 \ 
± y v 7m 2mc 2 J ' 

hence, these two particles will chase one another with different linear 
velocities: the difference, 

leads to a drift 

A/ = 6tt 



7m S 



rc 2 Mc 

of one particle with respect to the other per one revolution. Then in 
order to attain a drift of some 1 Angstrom per one revolution in a 
low orbit about the earth, the gyroscopes should spin with an angular 
velocity 

Q w 10 5 R- 2 sec' 1 , 

where R (cm) is the characteristic radius of inertia of the gyroscope. 
Some similar dragging (or quasi-dragging) effects, but in a quantum- 
mechanical description, were discussed in the monograph (Mitskievich 
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1969). Some other problems related to reference frames, gravitation 
and dragging in quantum physics, see in (Schmutzer 1975). 

It is worth finally noticing still another important manifestation 
of the dragging phenomenon (not merely in mechanics but essentially 
in the electromagnetic field theory) which is connected with results of 
DeWitt (1966) supporting the conclusion about a unity of gravotimag- 
netic and gravitoelectric effects in general relativity, which should be 
then both considered to belong to the dragging phenomenon. These 
results consist of two parts: 

(1) When a conductor is placed into a superposition of electric and 
gravitoelectric (Schwarzschild type) fields, or is considered (in the pres- 
ence of an electric field) in a translationally accelerated reference frame, 
inside such a conductor there vanishes not the electric field, but a cer- 
tain combination of that and of the gravitoelectric field, 

eE-mG = 0, (3.5.18) 

e and m being the charge and mass of electron (inside the conductor, a 
"free" electron gas has to be considered). Hence, in the conductor there 
must exist an electric field exactly corresponding to the acceleration of 
the monad congruence which describes the co-moving frame of this 
conductor. 

(2) When a superconductor is placed into a superposition of a mag- 
netic and gravitomagnetic fields (one may consider, e.g., the rotating 
Earth's field — say, the Lense-Thirring or Kerr one), or the supercon- 
ductor belongs to a rotating reference frame, inside this superconduc- 
tor there vanishes not the magnetic field, but its combination with the 
gravitomagnetic field, 

eB + 2muj = 0. (3.5.19) 

So it is usually said that in a rotating superconductor a corresponding 
magnetic field is generated. 

We propose to write in the case of a superconductor a general co- 
variant equation 

d{mu + eA) = (3.5.20) 

(this is a general covariant form of eqs. (3.5.18) and (3.5.19); cf. 
also (3.1.3)). One may characterize this effect as dragging of electric 
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and magnetic fields by the corresponding gravitational — or inertial - 
fields, cf. below (in section 4.3) kinematic "source" terms which appear 
in Maxwell's equations in non-inertial frames, and which are analogous 
to the Coriolis and centrifugal inertial forces terms in the equations of 
classical mechanics. 

It is clear that this new type of dragging makes it easier to measure 
gravitational fields, since it reduces gravitational measuring problems 
to the corresponding electromagnetic ones. Some proposals along these 
lines were made, e.g., by Papini (1966, 1969), although no general co- 
variant exact formulation was ever proposed (see also a consideration 
of the superconductivity theory in general relativity by Meier and Salie 
(1979) and Salie (1986)). 

It is also obvious that our approach makes it possible as well to con- 
sider non-stationary gravitational fields (and non-uniformly accelerated 
and rotating frames). Thus we predict that when a gravitational wave 
overrides a superconductor sample, corresponding electromagnetic os- 
cillations are to be observable inside it, the measuring problem being 
reduced to that of detecting electromagnetic oscillations inside a super- 
conductor. This is the case when no time-like Killing congruence could 
be present, which in the general opinion is however an indispensable 
premise for consideration of the dragging phenomenon itself. So, losing 
the symmetry in the sense of time, do we really lose the very thread of 
our argument?! 

It was however so easy to lose this tiny thread which does exist 
in highly idealised cases only. Adding merely a "grain" of impurity 
to the pure symmetry (here: stationarity), we extinguish the latter 
completely. But if a physical effect — or phenomenon — is so very 
vulnerable, it either does not exist at all, or its formulation is to be 
crucially revised. For example, in the Kerr field we encounter a great 
multitude of dragging effects; so, if from some distant source an in- 
finitesimally weak gravitational wave comes to this Kerr field region, 
could this lead to a total and instantaneous breakdown of all such ef- 
fects? Is it not much more likely that our definition of the dragging 
phenomenon (however natural it might seem to be) is in fact doomed 
to a radical revision? A symmetry is usually a good tool which makes 
it possible to dig up some principal results in the quickest and most 
elegant way, but we have not to reject a possibility of generalizing 
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these results beyond the limits of that symmetry approach, especially 
if the latter is so critically unstable. Description of the dragging phe- 
nomenon becomes very vague without the use of ideas of symmetry, 
because we have then no preferred reference frame — a very typical 
situation in general relativity. It is however quite plausible that the 
effects considered above in a stationary case, should hold also for fields 
essentially dependent on time. Thus in a cavity inside superconducting 
medium, there should be generated electromagnetic oscillations when 
a gravitational radiation pulse would override that region, so that the 
oscillations should imitate that pulse. Such a feature (one may speak 
on a tendency) is characteristic for gravitation itself (Mitskievich 1983). 
Somebody could insist that the notion of dragging is foreign to such 
effects, and probably it is not so important what name to give to the 
phenomenon the boundary of which is so diffuse. We are inclined nev- 
ertheless to speak about generalized manifestations of dragging, — let 
it be, for example, a combination of gravitoelectric and gravitomag- 
netic fields and corresponding effects. In fact, to every electromagnetic 
effect should correspond a certain gravitational effect (though not vice 
versa), which one may call a gravitoelectromagnetic effect, and a class 
of such effects pertains, in the limiting case of stationary space-times, 
to the dragging phenomenon area. 



Chapter 4 



The Maxwell field equations 



4.1 The four-dimensional Maxwell equa- 
tions 

To the end of the future applications of Maxwell's equations, we give 
here different representations of these equations in the four-dimensional 
form (without considering specific reference frames), with some relevant 
comments. 

We recall first the structure of the Lagrangian density (see Landau 
and Lifshitz 1971; Synge 1965; Mitskievich 1969) 



which yields Maxwell's equations through the standard variational prin- 
ciple. This Lagrangian density can be reshaped by singling out a di- 
vergence term: 



F^F"" = {2A V F^).„ - 2A v F' M, . li = *d(A A *F) — A ■ 5F. (4.1.2) 



The divergence (as usually) does not contribute to the Lagrange-Euler 
equations, but the variational principle should now be regarded as ap- 
plied to a Lagrange density which includes higher (second) order deriva- 
tives. The marvel of this reshaping is that the corresponding action 
integral can be varied according to the Palatini method known in the 
gravitation theory (Palatini 1919), i.e. when the both electromagnetic 




(4.1.1) 
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four-potential and field tensor are varied independently. Then the vari- 
ation with respect to the four-potential yields the Maxwell equations, 
and that with respect to the field tensor, the standard expression for 
this tensor as a four-dimensional curl of the four-potential (c/. the sit- 
uation in the gravitational field theory). As it is in the gravitational 
case, one has herewith to assume that Lagrangian densities of the other 
fields do not contain the electromagnetic field tensor, so that if it en- 
tered those Lagrangians before the reshaping (4.1.2) (which would be 
by itself exotic enough), then after passing to (4.1.2), one has to use for 
it in the Lagrangians four-curl of A a in its explicit form. Meanwhile 
the electromagnetic field Lagrangian does not explicitly contain in this 
formulation any derivatives of the four-potential. Thus the Palatini 
approach ascribes to the expression of the electromagnetic field tensor 
(2-form) 

F = dA (4.1.3) 

a dynamical sense. 

Now Maxwell's equations can be written in several equivalent forms: 

5F = 4tt7, dF = 0; 

F^, v = F^. v = 0; 

(V^gFn,, = -4*^=0?", (V^gFi?),* = 0; 
d* F = -47T *j, 5 * F = 0. 

Here the left-hand column represents variants of the Maxwell equations 
proper (let us denote all of them as (4.1.4a)), while the right-hand one, 
(4.1.4b), exspresses identical relations if the connection between F and 
A was imposed initially (F = dA), otherwise this is an independent 
system of equations. It is worth mentioning that asserting F to be a 
closed form, is not in general equivalent to consider it as an exact form; 
hence the theory is ramifying at this point. 

If one has assumed (4.1.3), only the equations (4.1.4a) have to be 
considered, these taking the form 

AA - d5A = 4irj (4.1.5) 

where it is convenient to impose the Lorenz gauge condition, 



(4.1.4) 



5A = 0. 



(4.1.6) 
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It is worth reminding that this condition is due to L. V. Lorenz of Copen- 
hagen, Denmark, and not to H.A. Lorentz of Leyden, Holland (see Pen- 
rose and Rindler 1983). Now the de-Rhamian A in the left-hand side 
includes the Ricci curvature tensor, see (1.2.49). 



4.2 The electromagnetic stress-energy ten- 
sor and its monad decomposition 

The Noether theorem (see Noether 1918, Mitskievich 1969, Mitskievich, 
Yefremov and Nesterov 1985) yields definition of the energy-momentum 
tensor density: 

X = yf^gT^ = -2- — . (4.2.1) 

We substitute here the Lagrangian density (4.1.1) using in it a tensor 
density of the weight |, 7 CTT = {—g) 1 ^g aT : 

„ 1 err 1 a/3 f f 



16tt 

An intermediate step will be 

1 di Tl3 
cm 4tt aT 0/37 d 9lxv 

where 

< .M/4 (_u!Lt-B\v 1 



Remark that g^v^ — = which results in vanishing of trace of the 
tensor T cm . Finally we come to the standard expression 

T cm u , = ~ (f, x F vX - \5;F KX F^ , (4.2.2) 

which can be written in a more symmetric form with help of the "crafty 
identities" (1.2.6): 

T cm ; = -^( F ^ FUX + F; x Ff). (4.2.3) 
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Due to (3.1.20), this expression is equivalent to 

T em ; = ~ (f» x F^ + F; X F^ = -i- F, X F» X (4.2.4) 

which immediately yields T cm ^ = 0. 

The following identities (cf. Wheeler 1962) hold for the electromag- 
netic stress-energy tensor: 

T C ^T C J U = (8n)- 2 [(B 2 - E 2 ) - A{E . B) 2 ]5 x m u, (4.2.5) 

or equivalently 

T cm ;T cm x u = (8n)- 2 [(B 2 + E 2 ) - A(E x B) 2 ]5 x m u, (4.2.6) 

In order to express the electromagnetic stress-energy tensor through 
the observables, let us consider 1-form T cra v ^r u dx^ (the electromagnetic 
stresses are to be considered separately). Then, making use of (3.1.9) 
and (3.1.14), we obtain 

T^dx^ = (8tt)- 1 (F mAj E a - F; x B x )dx» (4.2.7) 

Since 

F^V = *(£A*F) and F* x B x dx» = - * (B A F), (4.2.8) 
we have 

T^r^dx^ = (87T)" 1 * (E A *F + B A F). (4.2.9) 

After substituting here the decomposition (3.1.15) and its dual conju- 
gate (3.1.18) and easily identifying the scalar and vector products (see 
(2.2.15)), we obtain finally 

Tenfrda? = (8tt)- 1 [(^ 2 + B 2 )t + 2Ex B\. (4.2.10) 

It is clear that the component along the physical time of the reference 
frame (r), is the electromagnetic energy density, while the part lying 
in the three-space of the reference frame, the electromagnetic energy 
flux density (the Poynting vector), or the numerically identical with it 
(since c = 1) electromagnetic momentum density. We see that all these 
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quantities are expressed through the three-dimensional (but neverthe- 
less generally covariant) electromagnetic field observables (the electric 
field strength E and magnetic displacement B) in the very same way as 
they are expressed in the Maxwell theory in flat Minkowski space-time. 
This compels us to conclude that on the algebraic level, the electrody- 
namics theory has one and the same formulation in both special and 
general relativity, neither the gravitational field, nor the non-inertiality 
of a reference frame giving any essential contribution in its structure. 

As to the purely spatial part of the stress-energy tensor, we have to 
consider it not as an exterior form, but as a tensor, 

Stress = 6" <g> &„, (4.2.11) 

where b^ and b^ are semi-coordinated spatial bases (2.2.12), this tensor 
being (for convenience) of a mixed variance (once covariant and once 
contravariant). A combination of (4.2.11) and (4.2.3) yields 

Stress = (tory^F^V ® F vX b v + F; x b» ® F"X)- (4.2.12) 

In order to translate this to the reference frame language, we use a 
simple identity Something A = Something" 5*, with a subsequent inser- 
tion of 5 X = t x t k + b x , cf. (2.2.1). The terms with r yield the field 
vectors E and B via (3.1.9) and (3.1.14), while the terms with b x are 
expressed through the Levi-Civita axial tensor. We use here the ex- 
pression (3.1.11) as well as its natural counterpart which follows from 
the similarity of expressions (3.1.15) and (3.1.18). Further the second 
identity in (2.2.14) yields the final expression for the Stress: 

Stress = (Stt)- 1 ^ 2 + B 2 )W ®b ll -2(E®E + B® B)}, (4.2.13) 

the same as in the Minkowski space-time (cf. Synge 1965, p. 323). 

On the basis of invariants (3.1.16) and (3.1.17), the standard classi- 
fication of electromagnetic fields follows, consisting of the three types, 
electric (F^F^ < 0), magnetic (F^F^ > 0), and null (F^F^ = 0) 
ones, as it was the case for the Minkowski space-time. Like it was the 
case also in the flat world, such electromagnetic fields do not in general 
reduce to purely electric or purely magnetic fields (they are in general 
not pure fields) in no specific reference frame. However, also as in the 
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Minkowski world, the pure field condition is expressed as vanishing of 
the second invariant of electrodynamics, F* V F^ U = 0. When this con- 
dition is fulfilled, it becomes possible to find such reference frames in 
which the only non-zero observables are either the electric field strength 
(if F^ U F^ V < 0), or the magnetic displacement (if F^F^ > 0), or we 
are dealing with a pure radiation-type field (the both invariants are 
equal to zero), though this last case includes also fields which cannot 
be considered as propagating waves (e.g., some stationary fields; cf. 
(Synge 1974)). 

It is remarkable that in the last (pure radiation-type) case, an ana- 
logue of the Doppler effect holds: It is possible (in non-wave situa- 
tions too) to completely transform away the electromagnetic field by 
switching to moving reference frames, although this can be done only 
asymptotically (v — > c). 



4.3 Monad representation of Maxwell's 
equations 

First let us write down the action of differential operators on E, the 
frame-spatial orientation of the resulting quantities being conserved, if 
these quantities are non-scalar ones. The differential operations will 
be: The Lie derivative, 

£ T E = E a , x r x W + E x D x b° + E x A cX b°, (4.3.1) 

divergence, 

div E = 5E - G • E = E a . a -GmE, (4.3.2) 

and curl, 

curlE = *(t AdE). (4.3.3) 

These operations act similarly on the vector B. 

Consider now the Maxwell equations (4.1.4a), 5F = 4irj, taking 
into account (3.1.15): 



5F = S[(E A r) + *(B A r)] = 4tt/. 



(4.3.4) 
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Differentiation in the first term yields 

5(E A r) = E a . a r - E a , x r x dx <T + E a D a p dx p - D a a E + E a A aj3 dx^ 

= (div E)t - £ T E + 2E a D%dx f> - D^E, (4.3.5) 
and in the second term, 

5 * (B A r) = *d(B At) = curl B + G x B - 2[B • uj)t. (4.3.6) 

In the right-hand side of (4.3.4) we have to make use of the decom- 

(r) 

position of the four-current (3.2.8), j — u p(r + v), and to introduce 
the standard notations for charge and three-current densities, 

W W (r) 

p = u p, j =u pv, (4.3.7) 

which are non-invariant under transitions between different reference 
frames. We obtain then readily the Maxwell equations with sources, 
the scalar 

M 

div E = 4tt P +2uj • B (4.3.8) 

and the vector one, 

(3) 

curl B + Gx B= (£ T E - 2E l/ D^dx» + D a a E) + 4tt j . (4.3.9) 

The Maxwell equations without sources (4.1.4b), dF = (or, equiv- 
alently, 5 * F = 0), split relative to a reference frame in a complete 
analogy to the equations (4.1.4a), if we take into account that the ex- 
pression (3.1.19) differs from F (3.1.16) in (4.3.4) by an exchange of E 
by (-B) and B by E. Hence the scalar equation takes the form 

div B = —2uj • E, (4.3.10) 

and the vector one, 

curl E + G x E = -(£ T B - 2B v D v ll dxf t + D%B). (4.3.11) 

We see that in the presence of a gravitational field, and in a non- 
inertial reference frame characterized by acceleration, rotation, and 



80 



CHAPTER 4. THE MAXWELL FIELD EQUATIONS 



rate-of-strain tensor (thus also in the flat Minkowski world in such 
a reference frame), the Maxwell equations differ from the customary 
Maxwell's equations written for an inertial reference frame in the flat 
world, by the presence of additional terms. Such terms can be in- 
terpreted as supplementary effective charges (including an effective 
"magnetic monopole" charge), as well as the corresponding currents, 
cf. (Mitskievich 1969, Mitskievich and Kalev 1975, Mitskievich, Yefre- 
mov and Nesterov 1985). This is in a complete agreement with the 
known theoretical results claiming appearance of effective electric and 
magnetic (monopole) charges in rotating reference frames in the pres- 
ence of magnetic and electric fields correspondingly (see (Dehnen, Honl 
and Westpfahl 1961) and (Honl and Soergel-Fabricius 1961) for an 
approximate approach). Such charges and currents have of course a 
purely kinematic nature, though their objective existence is confirmed 
by the characteristic structure of the corresponding electric and mag- 
netic fields revealed by their action on charged particles when consid- 
ered in these reference frames. Obvious examples of this interpretation 
are presented in the following pages where we consider exact solutions 
of self-consistent systems of gravitational and electromagnetic fields to- 
gether with charged perfect fluid. 

It is however worth considering the analogy between classical me- 
chanics and electrodynamics in more general terms. In the former 
theory, the electromagnetic forces acting on charged particles, con- 
tain characteristics of the latters at most algebraically (in the three- 
dimensional description, electric field even enters the equations of mo- 
tion as an inhomogeneity) . In the field equations the analogous terms 
are the sources (mostly inhomogeneity, but the presence of a conduct- 
ing medium is revealed by a term proportional to the electric field, 
i.e. we encounter the same algebraic dependence, now of the field 
variable). This analogy between forces and sources can be traced more 
strictly and formally in the Lagrange- Euler form of equations of motion 
(for the both mechanics and field theory): they correspond to differ- 
entiation of the Lagrangian with respect to canonical coordinates (in 
the field case, the electromagnetic four-potential). Moreover, both the 
forces and the sources stem from one and the same terms in the inter- 
action Lagrangian density (which becomes the mechanical interaction 
Lagrangian when point-like sources are considered). In fact, this is inti- 
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mately related to an obvious generalization of the third Newtonian law 
of mechanics to include interaction of mechanical particles and fields. 
Now, in mechanics a non-inertial frame is revealed by appearance of 
the forces of inertia (the rotational case is usually considered: the cen- 
trifugal and Coriolis forces). One has also to expect an appearance 
of some (let us say) "sources of inertia", these however never having 
been considered in a far-reaching analogy with the forces of inertia of 
mechanics. But after we have come to the equations (4.3.8) — (4.3.11), 
this analogy became not a mere possibility, but an inevitable necessity. 

The Maxwell field equations seem to have a form much distinct from 
that of the equations of motion in mechanics, but it is now obvious that 
the terms connected with rotation of the reference frame (the vector a;) 
in (4.3.8) and (4.3.10), are direct analogs of the Coriolis force in me- 
chanics. Equations (4.3.9) and (4.3.11) contain the kinematic "sources 
of inertia" due to the translational acceleration and to deformations 
(dilatation and shear) of the reference frame; similar terms are known 
(though less discussed) in mechanics too. The problem now is to find 
out the role and importance of the corresponding non-inertial effects 
in the field theory, primarily in electrodynamics. The following pages 
contain some examples of how to handle this problem. 

4.4 A charged fluid without electric field 

We consider here exact solutions of the Einstein-Maxwell equations 
and equations for a charged perfect fluid, which were obtained as a 
generalization of the Godel space-time though having a more restricted 
isometry group (Mitskievich and Tsalakou 1991). These two families 
of solutions are described by a general metric 

ds 2 = e 2a (dt + fdxf - e 2f3 dx 2 - e^dy 2 - e 25 dz 2 (4.4.1) 

For the first of the families, (A), one has e a = a, e 2/3 = Eexp(2Bz + 
C) + Xz + v, 2 7 = 2Bz + C, 5 = (7 - (3) In a, / = bz/a 2 , while the 
invariant energy density of the fluid is 

p = ( /C a 2 )- 1 [(6 2 /a 2 - Ka 2 k 2 /2n) exp(-2Bz - C) - 3B 2 E], 
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its pressure p = B 2 E/(K,a 2 ), the invariant charge density 

p = -V2bk/(4ira 2 ) exp(-2Bz - C), 

where A = (kjW/tt - b 2 /a 2 )/2B, v = (A + D)/2B, B, C, D and 
E being integration constants and k a parameter which controls the 
switching on and off the charge and electromagnetic field. 

For the second family, (B), 2a = 2Bz + C, e 2/3 = (b/2B) 2 e - 2a - 
K(akz) 2 /2rr + 2Dz + E, e 2 ~< = F 2 , 5 = a-f3 + -f, f = -(b/2B) e - 2a , 
while the invariant energy density and pressure of the fluid are 

H = (2KF 2 y 1 e- 2a [2BD + Ka 2 k 2 (l - 2Bz)/2n], 

p = (2KF 2 y 1 e' 2a [2BD + Ka 2 k 2 (l - 2Bz)/2n], 
and the invariant charge density of the fluid 

p = -^[2abk{^F 2 Y x e- Za . 

Here B, C, D, E and F are integration constants. The both families of 
solutions contain also arbitrary constants a and b. For a special choice 
of the constants, the electromagnetic field is switched off, and the Godel 
(1949) cosmological model is recovered. 

We shall not discuss here the method of obtaining these solutions, 
but their hydrodynamical and electromagnetic properties are of con- 
siderable interest. It is natural to employ in their study the co-moving 
(with respect to the fluid) reference frame described by the monad field 

T = u = e a (dt + fdx). (4.4.2) 

Since the problem is a stationary one, the rate-of-strain tensor van- 
ishes for the both families of solutions, but acceleration and rotation in 
general survive, 

G = -a'dz, u = -f'dy. (4.4.3) 

Substituting into these expressions the functions entering the families 
(A) and (B), immediately shows for the family (A) that fluid is in a 
geodesic motion (G = 0), while its rotation is homogeneous in the 
whole of space (cu = {b/2a 2 )dy being an exact form), the rotation axis 
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directed along the y coordinate line. The family (A) fluid may be 
considered then as a superposition of two distinct fluids, a standard 
electrically neutral Godel fluid (if no cosmological term is introduced, 
it is a stiff matter), and a charged non-coherent dust. In contrast to 
this case, the fluid in the family (B) not only rotates (its rotation is 
inhomogeneous everywhere, although the rotation axis is still oriented 
along the y coordinate line), but moreover its particles are accelerated 
(non-geodesic motion). This is connected with non- homogeneity of the 
pressure (as opposite to the family (A) case) and by no means with 
any interaction with the electromagnetic field (c/. Section 3.2). The 
acceleration is directed along z axis, and it behaves exponentially with 
respect to this coordinate. 

For the both families (A) and (B), the electromagnetic field is de- 
scribed (in the coordinated basis) by one and the same potential 1-form 

A = ka(dt + V2z dx). (AAA) 

It is remarkable that this field is of magnetic type since 

F^F^ = 4a 2 k 2 exp(—2[3 — 25) > 0, (4.4.5) 

and moreover, it belongs to the pure type (F* u F^ y = 0); in the reference 
frame co-moving with the fluid, electric field is absent completely, in 
spite of the fact that the fluid is electrically charged (and the sign of the 
charge does not change anywhere). The magnetic field in the reference 
frame under consideration is 

B = V2ak ex-pin S)dy. (4.4.6) 

We have thus come to an apparently paradoxical situation in which a 
charged fluid does not generate (in a co-moving reference frame) any 
electric field at all, so that on its particles no lines of force of this field 
do neither begin nor end, although when these families of solutions were 
obtained, the Maxwell equations with a non- zero four-current density in 
the right-hand side were considered (and they were satisfied indeed); for 
some other examples see Islam (1985). It seems that this fluid generates 
a magnetic field only, and this does not influence the motion of the 
particles of this fluid, since the Lorentz force is identically equal to zero 
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(in the co- moving reference frame the particles are at rest by definition). 
The apparent paradox is resolved very simply (see Mitskievich and 
Tsalakou 1991) when it is taken into account that Maxwell's equations 
should be considered here in a rotating reference frame, taking the 
form (4.3.8) — (4.3.11). Then, a direct substitution of the quantities 
characterizing the fluid and electromagnetic field, into e.g. equation 
(4.3.8), shows that there occurs at all points an exact compensation 
of the charge density by the scalar product of the angular velocity of 
rotation of the reference frame and the magnetic displacement vector 
B, this justifying vanishing of the left-hand side of (4.3.8) which takes 
place simply by virtue of absence of the very electric field. Similar 
considerations hold concerning all other Maxwell's equations in this 
reference frame. 

In this problem, the crucial fact is that we consider an exact so- 
lution of a self-consistent system of equations for the fields and fluid, 
where the reference frame is objectively determined by the very state- 
ment of the problem, so that the final result asserting the possibility 
of a charged fluid which effectively generates no electric field, becomes 
obvious from the fact of no electromagnetic interaction between the 
(non-test) particles of this fluid. And the solution is an exact one! The 
particles of this fluid interact with each other only gravitationally (in 
the solution B, through pressure too). This effect has however a nature 
much simpler that the general relativistic one: it may arise in special 
relativistic and even non-relativistic circumstances. In such a highly 
non-linear theory as general relativity, it is difficult to immediately con- 
clude whether the fields (both gravitational and electromagnetic ones 
- though Maxwell's equations remain linear in general relativity too) 
are generated by the system of charges (and masses) under consider- 
ation, or there is a superposition of such a generated field with a free 
one. This is a case when general relativistic and reference frame con- 
siderations cannot help, thus it is worth treating this problem on the 
Minkowski flat space-time background. A readily obtainable conclusion 
is that in fact the both fields are present, that generated dynamically 
by the charged fluid, and a free (sourceless) field which is essentially 
magnetic and orthogonal to the three-velocity of the fluid particles. It 
is easy to see that even in the case of slowly moving fluid, i. e. when it 
is non-relativistic (though the electromagnetic field always represents a 
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relativistic object), such systems of charges and electromagnetic fields 
are realizable for which the Lorentz force completely vanishes. In the 
particular case of our generalization of Godel's solution, we would char- 
acterize it as a superposition of the charges Godel and Bonnor-Melvin 
universes. 



4.5 An Einstein-Maxwell field with kine- 
matic magnetic charges 

Let us consider now another example of non-inertial effects in exact 
solutions of Einstein-Maxwell equations. This will be a vacuum solution 
without any material carriers of charges (including magnetic charges, 
or monopoles), but in a rotating reference frame there will appear a 
distribution of purely kinematic monopoles, i.e. we shall have divB ^ 
0. 

To this end we consider the same form of metric as in the previous 
section, (4.4.1), but employ another Ansatz for the electromagnetic 
field, 

A — kt dy, (4.5.1) 

which automatically leads to satisfaction of Maxwell's equations with- 
out sources. Thus we may consider the case of absence of any matter 
in our new space-time. The electromagnetic invariants are now 

and \ (4.5.2) 

f; v f^ = o. J 

so that electric and magnetic fields should be mutually orthogonal. For 
simplicity it is natural to admit that the first invariant is also equal to 
zero, i.e. we would be dealing with a null electromagnetic field. Thus 
we postulate 

f = e 2 ^- a \ (4.5.3) 
The corresponding stress-energy tensor is 

Tem = ^ e -2(«-7)( (o) _ ( (o) _ 0(D). (4.5.4) 

4"7T 
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although a null tetrad is more appropriate in this case, we shall use the 
same orthonormal basis as in the previous section. 

Since \i = (no perfect fluid is present), we consider the case a + 
(3 = (this is obvious from the 00 and 11 components of Einstein's 
equations). A suitable choice of the coordinate z yields 7 = 5, and we 
come easily to the following solution: 

e 2a = e -2/3 = f -i = i^ z 2 + Az + Bjl = S = az + b. (4.5.5) 

Then the squared interval takes a remarkably simple form (let a = = 
b), 

ds 1 = e 2a dt 2 + 2dt dx - dy 2 - dz 2 , (4.5.6) 

with null x coordinate lines (or, which is here the same, null hyper- 
surfaces t = const), - - a special case of the pp- waves (see Kramer, 
Stephani, MacCallum and Herlt (1980), p. 233ff.). 
If we take r = = e a dt + e~ a dx, we get 

G = -ka'e-W 3 \ to = i/' e -^ (2) , D = 0, (4.5.6) 

as in the previous section, cf. (4.4.3), and 

E = -ke- a ~ 1 'd {2 \ B = -ke-^tf^ . (4.5.7) 

Turning now to Maxwell's equations, we find that u> • B = 0, so 
that divE = 0, but u • E = \kfe~ a - 2 ~i ^ 0, hence divB ^ too, and 
we come to a distribution of effective magnetic charges as kinematic 
manifestation of the reference frame rotation. The same conclusions 
immediately follow from the observations concerning orientation and 
^-dependence of the (co)vectors E and B. As to the observability of 
such purely classical kinematic monopoles, we insist that it is perfectly 
possible, although there is still some deficiency in the definition of the 
physical three-space of a rotating reference frame, since this submani- 
fold is only local (non-holonomic) , so that it is impossible to speak un- 
equivocally about integration over its parts (in order to produce electric 
and magnetic fluxes). But if infinitesimal regions are considered (and 
for a continuous distribution of magnetic charges this is just the case), 



4.5. KINEMATIC MAGNETIC CHARGES 



87 



it is quite plausible that such kinematic charges have to be as real as, 
e.g., the Coriolis force is. 

Thus we conclude that when experiments on detection of dynami- 
cal magnetic monopoles are performed, it is important to clearly dis- 
tinguish them from kinematic magnetic charges. However when decay 
processes are considered, no such a problem could arise, though at a 
classical level influences of the both seem to be largely identical (in the 
kinematic case, the very rotation of the Earth may simulate a presence 
of monopoles). 
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Chapter 5 

The Einstein field equations 



5.1 The four- dimensional Einstein equa- 
tions 

Let us begin with some remarks on the structure of four-dimensional 
(priory to performing the monad splitting) Einstein's equations. In 
chapter 4 we mentioned the Palatini method when the Maxwell equa- 
tions were deduced using the electromagnetic Lagrangian. In the grav- 
itation theory, due to the principle of equivalence, there exists only a 
relative gravitational force acting on test particles, to which a relative 
field strength corresponds (we would rather prefer to use the expression 
"gravitational inhomogeneity field", since observable quantities are in 
any case determined relative to a reference frame). In order to come to 
this concept, one may use the geodesic deviation equation, 



which is deduced from the assumption of geodesic motion of free test 
particles, V M w = 0, and includes a vector field w commuting with 
the four-velocity field u of the particles, so that w describes closeness 
of world lines in such a time-like congruence (hence, this equation is 
related to the focusing effects of the motion of particles, and to the 
concept of space-time singularities). From (5.1.1) it is clear that the 
role of the gravitational inhomogeneity field is played by the Riemann- 
Christoffel curvature tensor (this is the only ingredient in the geodesic 
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deviation equation which does not depend on characteristics of mo- 
tion of the particles themselves). If one applies the electromagnetic 
Lagrangian (4.1.1) while throwing off a divergence term according to 
(4.1.2), it is natural to find an analogy between the quantities 

A g a/3 , 5F^R\ p<J = R aP . (5.1.2) 

Thus we come to the gravitational Lagrangian density 

Z g = ^g^R aP (5.1.3) 

in full conformity to the electromagnetic theory, while the constant 
coefficient involves the Newtonian gravitational constant 7 (the Ein- 
steinian constant is also widely used, k = 8777). Every variational 
method, including that of Palatini (1919) (c/. for other approaches 
(Landau and Lifshitz 1971, Fock 1964) and (DeWitt 1965) where an 
important concept of the second variational (functional) derivative is 
also considered), now yields the well known Einstein gravitational field 
equations, 

R a p - ^g a /3R = -nT a p. (5.1.4) 

Hence it is logical to expect that the monad splitting of these equations 
(supplemented also by the sourceless equations in analogy with the 
Maxwell theory), should yield scalar and vector field equations which 
would exhibit some analogy with Maxwell's equations under the simi- 
lar splitting (there arise also three-tensor gravitational field equations 
which are of course sui generis). 

Any physical field has certainly its own characteristic features (oth- 
erwise is would not deserve the status of a separate field); this is true 
with respect to the gravitational field too whose equations have their 
own peculiarities. We shall disclose the latters in course of the following 
consideration, but the problem is so extensive and complicated that it 
is worth being a subject of separate study in another publication: here 
we are interested mainly in applications of the monad formalism. 
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5.2 Monad representation of Einstein's 
equations 

First we consider equations containing no source terms which are similar 
to Maxwell's equations (4.3.11) and (4.3.12). These turned to be iden- 
tities when the electric field strength and magnetic displacement were 
expressed through the four-potential (the reader may do this easily in 
the monad language as an exercise); the new gravitational equations 
are identities in the same sense. Remember expressions for divergence 
and curl of section 2.4, in particular (2.4.4). Inserting them into (2.4.7) 
and (2.4.9) (the first of these relations is repeatedly used in the second 
one), we obtain 

divu = G»u (5.2.1) 

and 

curl G = 2 (£ r w - 2u a D al3 dx 13 + Quj) . (5.2.2) 

When these equations are compared with Maxwell's equations (4.3.10) 
and (4.3.11), one sees that they are in fact the same up to a substitution 
of G instead of E and (— 2uj) instead of B; this substitution is exactly 
the same which can be observed in (3.3.11) and the term analogous to 
G x E of (4.3.11) naturally vanishes. 

As to Einstein's equations proper, they enter the monad construc- 
tions via Rap^r 13 , RapT^ 13 , R a pb a <g) y 3 , while it is the best to denote 
T a /3T a T^ =: e, Tap^b 13 =: a, T a/3 b a (g> b 13 =: rj, where e is energy den- 
sity scalar, a energy flow density (or equivalently, momentum density) 
covector, and rj tensor of the purely spatial stresses (or, momentum 
flow density), all taken with respect to the reference frame r (so that 
for perfect fluid in a co-moving frame, r = u, one has e = ji, a = 0, 
V — p(9 ~ u ® u ), c f- (3.2.6); for perfect fluid in an arbitrary frame, e, 
<j, and rj can be read off from (3.2.7)). Then energy-momentum tensor 
takes the form 

T = er®T + T®a + a®T + rj ) (5.2.3) 

while r ■ a = 0, so that the contracted energy- momentum tensor is 
trT = e + tr?7. 
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The obvious identity 

i^rV = - r a mf3 y 

can be easily brought into the form 

i^r V = 2Q, a r a - G a , a + D al3 D a/3 - A aj3 A al3 , (5.2.4) 

and it yields finally a scalar equation (its special case for G = is the 
Raychaudhuri equation, see (Ryan and Shepley 1975)): 

div G = ^(e - trrj) + 26, a T a - G«G-2cj«cj + D a/3 D a ^. (5.2.5) 

Here we used definition of divergence, (2.4.5), and the fact that A a pA alS = 
2cj • a;. If now the relation G ,a A Q/3 (ix /3 = w x G is inserted into the gen- 
eralized Codazzi equations (2.4.12), one obtains the vector equation 

curl uj = -kg + 2uj x G + (D£. v l% - 2Q >a )b a . (5.2.6) 

The tensor equation can be obtained without difficulty on the basis 
of the generalized Gauss equations (2.4.13), using as well relation 

Ra/3 = R K al3\T X TK + R K a/3xb^ 

where the left-hand side (Ricci tensor) is to be expressed through the 
energy-momentum tensor (sources of Einstein's equations), while 

so that this purely three-spatial quantity is easily expressible through 
the reference frame characteristics only: 

Ra^r 5 = l -G^(b^ + b^ a ) - £ T D, a - G a G a + D^D v a 

+u a u^ + u • uj b ay + 2e flS ( a D s l) uj' 1 . 
At the same time, the crafty identities (1.2.9) yield 

R a( 3 7 sK b x b l b i = -Rc t p 1 5T l3 T 5 e KX a e ilu 1 - ^R{b^b Xv - b KU b Xfl ) 
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-R ay e KXa e^ + RpT^Tsib^bxv - b KI/ b XfJ ), 

but this is exactly the quantity which should be substituted into the 
right-hand part of the Gauss equations (2.4.13), or in their contracted 
form (2.4.14), which makes no difference whatever since for a three- 
dimensional case the Riemann-Christoffel and Ricci tensors are equiv- 
alent (they both possess one and the same number of independent non- 
trivial components). The subsequent calculations represent a mere rou- 
tine, yielding finally 

r AM = \G K .,Mbl + - £ T D Xfl - G X G, + 2D X „D; - 2QD Xfl 

+2e mx D 5 ^ + ^eb Xfl - K ( VXfl - -b^trrj). (5.2.7) 

This tensor equation can be found in literature (however, in forms dif- 
fering much one from the other, not only concerning notations; cf. 
(Vladimirov 1982, Zel'manov and Agakov 1989)). It has no counterpart 
in Maxwell's electrodynamics indeed, its characteristic feature being 
presence of Ricci tensor of the three-dimensional space of the reference 
frame, so that (5.2.7) may be interpreted also as means for obtain- 
ing this latter quantity. We shall not discuss here the important role 
of obtained equations from the point of view of the Cauchy problem 
approach. 

In the dynamical equations which include sources of the field, spe- 
cific character of gravitation manifests itself more clearly, although 
there are certain common terms existing in electromagnetic equations 
too. The most striking feature of (5.2.6) is absence of time deriva- 
tive of the acceleration vector (cf. £ T E in (4.3.9)). Dehnen (1962) 
succeeded in eliminating this distinction but at the cost of a highly 
artificial approach, and a coordinate representation was essential (cf. 
also our approach (Mitskievich 1969) which is probably more transpar- 
ent). One knows however from the formulation of the Cauchy problem 
for Einstein's gravitational field that the noted distinction of Einstein's 
and Maxwell's theories is a crucial one (see e.g. (Mitskievich, Yefremov 
and Nesterov 1985)). 

Thus a question arises, how far-reaching can be the natural anal- 
ogy between gravitation and electromagnetism (without abandoning 
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the standard context of Einstein's theory of gravity)? The study of 
the equations of motion of electrically charged test particles in chap- 
ter 3 has already shown that this analogy is very profound, starting 
with the four- dimensional form of these equations, (3.1.4), and being 
revealed more comprehensively in the monad splitted equations (3.3.9) 
- (3.3.11). In this connection there was even extinguished the basic 
distinction of one theory from another, the uniqueness of the gravita- 
tional equivalence principle, although the latter is manifested primarily 
by equations of motion of the particles. In equation (3.1.4) this enters 
as the repeated participation of the four-velocity (or the absence of a 
similar repetition for the four-potential). Equations of the very gravita- 
tional field are however much more involved than one would expect: if 
the set of equations, being in fact identities when acceleration, rotation 
and rate-of-strain are expressed through the monad field, still remains 
quasi-Maxwellian in its form, the more profound dynamical equations 
of gravity, (5.2.4) and (5.2.5), differ from the corresponding Maxwell 
equations, (4.3.8) and (4.3.9), radically enough. 

Does there exist any roundabout way which could outflank this 
difference without revising the already existing theory (for such a re- 
vision there exist no reasons at all according to the profound belief 
of the most of the specialists)? This way really does exist, but it 
is not as straightforward one as the approach just considered above. 
This alternative representation of the gravitation theory is based on 
the Riemann-Christoffel curvature tensor instead of that of Ricci used 
in Einstein's equations; moreover, field equations now are first-order 
differential ones with respect to this tensor (third-order to the met- 
ric tensor), while the Ricci tensor enters Einstein's equations purely 
algebraically. Nevertheless, the physical contents of this higher order 
theory is the same as it was for Einstein's theory. The gravitational 
field equations take then the form of contracted Bianchi identities with 
a substitution of the right-hand side of Einstein's equations instead 
of the Ricci tensor; then the four-dimensional equations (without in- 
troduction of any reference frame approach) show analogy with the 
Maxwell equations: 
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(old quasi-Maxwellian gravitational equations; the new ones are (5.4.3)). 
A detailed study of these equations see in (Mitskievich 1969) where also 
the variational formalism (with the corresponding Lagrangian density) 
was proposed. The first attempt of splitting the gravitational inhomo- 
geneity field Ra^s in a reference frame sense (not exactly the monad 
one) was made by Matte (1953) and discussed in a modernized form by 
Zakharov (1972) (see also Zel'manov and Agakov (1989)). Hereto the 
notations were introduced, 

X a( s = R^pT^r", Y a(3 = -R^pT^r" ', Z a(3 = -R* au pT^T u , (5.2.9) 

which are clearly similar to (3.1.9) and (3.1.15) for the electromagnetic 
case (though they are quadratic with respect to projections onto r). 
Dashes were written here since we shall use below somewhat different 
quantities X and Y: 

X a p = CftavpT^T" ', Y a/3 = —C * pT^T u . (5.2.10) 



5.3 The geodesic deviation equation and 
a new level of analogy between grav- 
itation and electromagnet ism 

As it was pointed out in section 5.1, the existing exact analogy between 
gravitation and electromagnetism enables one to formulate the gravita- 
tional Lagrangian along the same lines as it was for the electromagnetic 
one. The directing principle was the use of equations of motion of a 
particle interacting with the corresponding field in order to determine 
the field strength concept; for gravitational field one has to consider the 
geodesic deviation equation (5.1.1). Obviously, it is worth making use 
of this equation also for finding hints of how to introduce analogues of 
the electric and magnetic field vectors connected with the gravitational 
inhomogeneity field. 

Let us express the full Riemann-Christoffel tensor in (5.1.1) through 
the Weyl conformal curvature tensor, Ricci tensor and the scalar cur- 
vature. Consider now the term C a p 1 su l3 u 1 w s which may be written 
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also as C a/ 3 7 sdx al3 dx' yS (-,u,u,w). The Weyl tensor may be considered 
as the representative of a free gravitational field, the other terms en- 
tering the Riemann-Christoffel curvature being expressed through the 
gravitational field sources (cf. the ideas of Pirani and Schild (1961) and 
Wheeler (1993)). The quantities X and Y which enter the Weyl tensor, 
are not only symmetric in their indices, but also traceless, so that they 
are the only irreducible representatives of free (intrinsic) gravitational 
inhomogeneity field. We apply here the expression (2.2.21) for the basis 
2-forms and note that the necessary contractions are 



(r A &")(•, u) = \ u (TV^-b"), 

Z 

*(r A &")(•, u) = -\ u v x If, 
z 

1 M 

*(r A b^)(u, w) — —- u (v x wy. 
z 

Recall now that in the electromagnetic case one has 

1 (t) 



F(.,u) 



u 



[{E»v)t + E + {v x B)} 



(5.3.1) 

(5.3.2) 
(5.3.3) 

(5.3.4) 



which gives a combination of the Lorentz force with the work of the 
field pro unit time (and unit charge of the particle); the absence of the 
second vector, w, is here but natural. 
In the case of the Weyl tensor we have 



C a f3 y sdx a ^dx l5 (-, u, u, w) = 
where 



[(Af • v)t + M + (v x iV)] (5.3.5) 



M := 



X, 



a(3 



W V h 



and 



N 



+ -Y aP {v x wf 



(r) 



-X af3 (v x w) p + -Y aj3 [w v'-w 



b a 



b a . 



(5.3.6) 



(5.3.7) 



For X and Y see (5.2.10); it is clear that the analogue of Z but without a 
dash would merely coincide with X, so it is superfluous. The formulae 
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(5.3.4) and (5.3.5) are strikingly similar to each other; the factor | 
in (5.3.4) is due to the definition of the 2-form F, and the second 

(t) 

power of the standard relativistic factor u results from the repeated 
appearance of the four- velocity u in the left-hand side of (5.3.5). It is 
however impossible to completely identify M and N as exact analogues 
of the electromagnetic quantities E and B since the formers include 
a dependence on u, the velocity of particle, and on another vector 
field w (density of world lines of a cloud of test particles), so they 
cannot be interpreted as pure field quantities. Consider instead an 
equation generalizing the geodesic deviation equation (5.1.1) to the 
case of electrically charged test masses (with equal charge/mass ratio 
e/m), 

VuVuW = —(F af3 u^). n w^dx a + R(u, w)u. (5.3.8) 



We see that there exists certain difference in the structure of the electro- 
magnetic and gravitational (curvature) terms: the latter is quadratic in 
u. When expressed with respect to a reference frame, equation (5.3.8) 
contains alongside with the term independent of the three- velocity v 
(quasi-electric vector) and that linear in the velocity (in electrodynam- 
ics, the vector product containing the magnetic displacement vector B), 
also a term quadratic in the three- velocity of the particle, which cannot 
have any counterpart in the electromagnetic theory. As to the other 
vector field, w, it is specific for the deviation type equations, and it is 
present to the same extent also in the electromagnetic part of (5.3.8). 

The only reasonable way out of this situation seems to be an in- 
troduction of the already well known symmetric three-tensors, X and 
Y, present in the motion-of-the-particle dependent quantities M and 
N, (5.3.6) and (5.3.7). The first is clearly most closely related to the 
electric type and the other to the magnetic type field, but the intrin- 
sic gravitational inhomogeneity field tensor (the Weyl tensor) exhibits 
certain peculiarities when entering the geodesic deviation equation. It 
remains to interprete these as manifestations of the higher tensor rank 
of the gravitational field. 



98 



CHAPTER 5. THE EINSTEIN FIELD EQUATIONS 



5.4 New quasi-Maxwellian equations of 
the gravitational field 

Once contracted Bianchi identities (1.2.23) read 

R"L»;X = (5.4.1) 

or equivalently, 

R a p 1 5-,a = R/3 T ,s ~ Rps-a- (5.4.2) 

When the right-hand side of Einstein's equations is substituted into the 
right-hand part of (5.4.2), this yields the old quasi-Maxwellian equa- 
tions (5.2.8). It is however better to expess the Riemann-Christoffcl 
tensor in the left-hand side of (5.4.2) (or (5.2.8)) in terms of the Weyl 
conformal curvature, Ricci tensor and scalar curvature (see (1.2.24)), 
and further express the latters in terms of the stress-energy tensor and 
its trace using Einstein's equations. 

It is easy to find that (5.4.2) and (5.2.8) yield in this way new 
equations 

k 1 

C aT ^ u ' = — \T T v;ix T TjJj . l/ ~^{T^g T y T ]i/Q T ^)] (5.4.3) 

(in fact, one can find these equations already in the book by Eisenhart 
(1926), though without their physical interpretation). They hold for 
arbitrary solutions of Einstein's equations; moreover (see Lichnerowicz 
(I960)), if a solution of Einstein's equations is taken as initial bound- 
ary values, it is reproduced in the whole space-time as a solution of 
equations (5.2.8) (so of course also of (5.4.3)). These equations de- 
termine (when considered in analogy with Maxwell's equations for the 
field tensor F) the intrinsic gravitational inhomogeneity field tensor 
C (the other parts of the Riemann-Christoffel tensor describe another 
- extrinsic — kind of the gravitational inhomogeneity field which is 
algebraically reducible to its sources, the stress-energy tensor T). 

It is important that non-contracted Bianchi identities yield exactly 
the same equations (the contraction does not lead to any information 
loss in this case): they can be written with help of the Weyl tensor, 
as it already was the case for the contracted identities (c/. (5.4.2) and 
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then (5.4.3)). Thus the dual conjugated Weyl tensor now enters the 
left-hand side as C* Tflu ,a . But C* Tflu = C aT * u , so that the resulting 
quasi-Maxwellian equations are nothing else than the equation (5.4.3) 
with dual conjugation in the indices ji and v\ As a matter of fact, this 
is equivalent to representation of the Bianchi identities as vanishing of 
divergence of the right-hand side of (1.2.9). 

Now we have to express the left-hand side of (5.4.3) through the 
quasi-electric and quasi-magnetic space-like symmetric tensors X and 
Y . To this end we first substitute relation (2.2.21) into the coordinates- 
free expression of Weyl's tensor and then take account of the definitions 
(5.2.10). The resulting decomposition is 

^C aMS dx af5 <g> dx~< 5 = Xp S [(r A b 13 ) <g> (r A b 5 ) - *(r A iP) <g> *(r A b 5 )] 

+Yp 5 [(r A b 13 ) <g> (r A b 5 ) + *(r A b ) <g> *(r A b 5 )}. (5.4.4) 

An analogy with the definition (3.1.15) is here obvious, although sum- 
mations over the indices of X and Y components are shared among 
the both two-forms of the tensor (double bivector) basis. The pure 
component form of the same expression is 

Caf3~ ( 5 = XpsTaT^ + X al TpT 5 — X^TpT^ — X^T^s — X^ E ^^E^xT^ T V 

+Y px E^ x r a T u - Y aX E l5 v %T v + Y M E aP "\T v - Y^E a/3 " x r s r u . (5.4.5) 

Consider now the left-hand side of the quasi-Maxwellian equations 
(5.4.3), i.e. four-dimensional divergence of the Weyl curvature. A 
somewhat lengthy calculation yields the following decomposition with 
respect to the monad r: 

C a/3 y8 ,a = -Tp(TyI s - r 5 I 7 ) - Tpll x e y s\ + Tulips - Tslllp-f + Wp X e~ f5x 

(5.4.6) 

with 

IIL3 7 e^ A - II A = 0, lll a pg^ = 0, \ 
Wp 7 e^ x + I A = 0, W a pg a P = J [b <} 

(consequences of the four-dimensional Ricci identities for C as well as 
of its tracelessness) . Here 

Is = X aa ^ 5 + G a X a5 + 3uj a Y aS + D av Y x e 5vX , (5.4.8) 
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Us = Y a ^b% + G a Y a5 - 3u a X a5 - D w X^e SvX , (5.4.9) 
5 1 

IH(/35) = £ T Xf)S — -(XpD aS + XgD a p) + -^(XpCaSy + X^e Q/37 ) 

+4e^+X^Z^ 5 -G'% A e feA -^ 

(5.4.10) 

IV m = £ T Y P s - \{Y$D a5 + Y s a D a p) + \j~\Y«e a&1 + Y s a e a ^) 

+46^ 5 +VD^^-G'^X A e^ A -G'^X 5 A e ^ A +ix KA;a (e a ^6 A5 +e a5K 6 A/3 ) 

(5.4.11) 

(it is clear from (5.4.7) that only the symmetric parts of III and IV are 
independent constructions). Since the both quasi-electric and quasi- 
magnetic gravitational inhomogeneity fields are two-index tensors and 
not vectors, as it was the case for electromagnetic field, it is inappro- 
priate to introduce here operations as curl or div (c/. Section 4.3), as 
well as to make use of three-dimensional scalar and vector products. 

In a vacuum, the field equations of gravitational field (the higher- 
order counterpart of Einstein's equations) reduce to vanishing of the 
quantities (5.4.8) - (5.4.11). A straightforward comparison with the 
corresponding parts of Maxwell's equations (4.3.8) - (4.3.12) shows 
existence of every parallels to electrodynamics, with natural differences 
in the coefficients, while the new specific terms in the gravitational 
field equations are those which cannot in principle arise when three- 
dimensional representation of fields is realized with vectors (E, B) and 
not tensors (X, Y). The only exclusion is appearance of terms with u 1 , 
the reference frame rotation vector, in (5.4.10) and (5.4.11), but this 
can be easily understood since this comes from the expression of the 
time-derivative, e.g. 

£ T Xj3S = Xps. a T a + X a $T a .p + Xp a T a .g, 

which contains as compared with £ T Ep, an extra term with r a -s, and 
the complete compensation (as it was in (4.3.5) and (4.3.1)) is now 
impossible. In order to facilitate such a comparison, we give here ex- 
pressions for Maxwell's equations of Chapter 4 in the component form: 



E a , a + G a E a + 2u a B a = 4tt p , 



(5.4.12) 
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B a , a + G a B a - 2u a E a = 0, (5.4.13) 

(3) 

£ T E S - 2E a D a5 + 2BE 5 + G»B x e 5uX + B^ a e aSK = -4tt j 5 , (5.4.14) 

£ T B 5 - 2B a D a5 + 2QB 5 - G"E x e SvX - E K > a e aSK = 0. (5.4.15) 

It is obvious that this representation of the gravitational field equations 
is completely in line with the form of Maxwell's equations, in contrast 
with the direct decomposition of Einstein's equations: in particular, the 
time-derivatives always enter together with the analogues of curls of the 
counterpart fields (this was exactly the stumbling-block in section 5.2). 

As to the right-hand side of the quasi-Maxwellian equations, we 
shall restrict our consideration to a perfect fluid in a co-moving ref- 
erence frame, with the energy-momentum tensor (3.2.6). Then the 
quasi-Maxwellian equations (5.4.3) take form 



K 



(5.4.16) 



k(h + p)uj 5 , 
111(0*) = n(fi + p) 



(5.4.17) 
(5.4.18) 

TV m = 0, (5.4.19) 

where for the left-hand sides, expressions (5.4.8) - (5.4.11) should be 
substituted. The right-hand side expressions are now remarkably sim- 
ple, and it is easy to check that they satisfy conditions (5.4.7) (since in 
(5.4.16) - (5.4.19) equations of motion (3.2.10) and (3.2.11) were taken 
into account, they should be used in this check also). 



5.5 Remarks on classification of intrinsic 
gravitational fields 

Now it is clear that gravitational inhomogeneity fields can be classified 
in analogy to the very simple classification of electromagnetic fields 
which divides the latters into the fields of electric, magnetic and null 
types (see e.g. Synge (1956)). Very naturally, this classification is 
completely in line with the well known algebraic classification of Petrov 
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(1966) (see also (Geheniau 1957, Pirani 1957, Bel 1959, Debever 1959, 
Penrose 1960; Kramer, Stephani, MacCallum and Herlt 1980)). 

There exist only four independent invariants of the Weyl tensor C: 



J-l — ^afiyS^ , 1-2 — ^afiyS^ ) 



(5.5.1) 



(they can be found as the identically non- vanishing constructions among 
the fourteen invariants of the Riemann-Christoffel tensor when the lat- 
ter is substituted by the Weyl conformal curvature tensor, see e.g. (Mit- 
skievich, Yefremov and Nesterov 1985), p. 76). When one expresses 
the Weyl tensor via (5.4.4) or (5.4.5), a representation of the invariants 
in terms of the quasi-electric and quasi-magnetic gravitational fields 
emerges: 



h = KXpsX? 6 - YpsY? 6 ), I 3 = 16X$(X$X2 - SY^YJ), 
h = -ISXpsYf*, h = 16Y£(YPY2 - 3X^X2) 

(5.5.2) 

Note that these invariants do not change under the both independent 
infinite-parameter groups of transformations, those of systems of co- 
ordinates and of reference frames. They are closely connected with 
(anti-)self-dual fields (cf. the electromagnetic case, section 3.1). To 



this end observe that 
where 
Then 

and finally 



^a(3-yS = =^ Cal3y5 (5.5.3) 



Ca^s'-— C al 3.ys T iC*^. (5.5.4) 



X a p--=Cani3v ^T V = X a/3 ± lY afi (5.5.5) 



h T ih =C a ^ s C lS tv C erj a/3 = 16 XPXJX* . (5.5.6) 



Similarly, 



I 1T U 2 = 8XPX% . (5.5.7) 

Since the quadratic constructions formed exclusively of the real (not 
(anti-) self-dual) X or Y are positive definite, it is clear that the sign 
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of Ii determines the electric (plus) or magnetic (minus) type of the 
gravitational inhomogeneity field, while I\ — corresponds to the 
semi-null type field. The gravitational inhomogeneity field is purely 
quasi-electric, purely quasi-magnetic, or null when all the other invari- 
ants vanish (in the case of electromagnetic field the situation is more 
simple since there exist only two corresponding invariants, and vanish- 
ing of the only one second invariant provides the condition of a pure 
field). When we conclude that the field is purely (quasi-)electric (or 
magnetic), this does not mean that the alternative field is absent iden- 
tically: it simply can be transformed away by a choice of the reference 
frame. The gravitational null case has the same propery as its electro- 
magnetic counterpart: the null field merely changes its amplitude under 
transitions between different reference frames (a manifestation of the 
Doppler effect in the non-wave representation), with a possibility to 
transform away (at least, locally) the whole of the field in an asymp- 
totical sense (when the velocity of a new reference frame is tending to 
the speed of light with respect to the initial frame). 

This amazingly complete and simple analogy between electromag- 
netism and gravitation is well known indeed, but its existence is worth 
being mentioned since it is widely underestimated. We shall not go 
into more detail concerning the Petrov classification in connection with 
the properties of the Weyl tensor invariants, as well as the correspond- 
ing properties of X and Y: these aspects may be considered as not so 
closely related to the reference frames problems. 



5.6 Example of the Taub-NUT field 



We consider here first the Taub-NUT (in fact, especially NUT) field 
since it incorporates the both quasi-electric and quasi-magnetic proper- 
ties which are already very well known (for references see e.g. (Kramer, 
Stephani, MacCallum and Herlt 1980)). Its metric reads 



ds 2 = e 2a (dt + 2/cos^0) 2 - e~ 2a dr 2 - (r 2 + l 2 )(d§ 2 + sin 2 ^0 2 ), 




(5.6.1) 
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Here e 2a = (r 2 —2mr—l 2 )/(r 2 +l 2 ). The only non- vanishing independent 
curvature components in this basis are 

fl(0)(l)(0)(l) = -#(2)(3)(2)(3) = 217, 

#(0)(2)(0)(2) = #(0)(3)(0)(3) = --R(l)(2)(l)(2) = --R(l)(3)(l)(3) = f 7 , 
#(0)(1)(2)(3) = -2i2( )(3)(l)(2) = -2i2( )(2)(3)(l) 

where 

U = (Z 4 + 3m/ 2 r - 3Z 2 r 2 - mr 3 )(r 2 + Z 2 )" 3 , 
V = l(ml 2 - 3Z 2 r - 3mr 2 + r 3 )(r 2 + I 2 )- 2 . 

The most convenient choice of the monad field is 

t = 6 ( ® = e a {dt + 21 cos M<P) 

which describes an accelerated and rotating (but without deformations) 
reference frame, 

G = -a'e a 6 {1 \ u = -le a (r i + l 2 )- 1 0W (5.6.4) 

(the both vectors have radial orientation, and an effective "spherical 
symmetry" takes place). Thus from the point of view of section 5.2, 
we can treat the NUT field as possessing analogues of the both electric 
(G) and magnetic (u) fields, their radial orientation corresponding to 
a similarity to the magnetic monopole field superimposed on a point 
charge field. 

The latter interpretation depends however crucially on choice of the 
reference frame, since it is always possible to consider a non-rotating 
frame, e.g. with 

r = N{dt + 2ldxf>), (5.6.5) 

N 2 = (r 2 - 2mr - l 2 )(r 2 + Z 2 )[(r 2 + 1 2 ) 2 - 4l 2 (r 2 - 2mr - I 2 ) tan 2 (0/2)], 
thus showing analogy with an electric field only. We do not present here 
the corresponding acceleration vector which has components along axes 
r and ■&, since it is ruther unsightly. 



2V, 



(5.6.2) 



(5.6.3) 
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In the rotating reference frame (5.6.3), the tensors of quasi-electric 
(gravitoelectric) and quasi-magnetic (gravitomagnetic) gravitational in- 
homogeneity fields (5.2.9) are diagonal with the components 




^(i)(i) = -2U, X {2){2) = X (3)(3) 

y ( i) ( i) = 2V, y (2)(2 ) = y ( 3) ( 3) = 

Then the Weyl tensor invariants read 

I l= 48(U 2 -V 2 ), I 2 = 9QUV, 

h = 96U(U 2 - 3V 2 ), I 4 = -96V(V 2 - 3U 2 ). 



(5.6.6) 



(5.6.7) 



Here the above remarks on the classification of gravitational inhomo- 
geneity fields (section 5.5) are clearly applicable. 

For the (anti-)self-dual Weyl tensor (5.5.4) and corresponding space- 
like complex tensors (5.5.5), it is easy to find that 

U±iV= -(r 2 + l 2 )- 3 (m t H)(r ± U) 3 (5.6.8) 

(cf. a superposition of electric point charge and magnetic monopole 
fields in Maxwell's theory). Then 



± 

X- 



( -2 

1 | (UTiV). (5.6.9) 

V o o 



It is worth mentioning that this type of solution may be obtained in a 
purely linear theory using cohomology and twistor methods, cf. (Hugh- 
ston 1979). 

As to the description of quasi-electric and quasi-magnetic gravita- 
tional inhomogeneity fields in a non-rotating reference frame, e.g. that 
of (5.6.5), it involves them both in the Taub-NUT space-time, although 
they enter not so symmetric as it is the case for the rotating frame (see 
(5.6.6)). The transition between the two frames is given by 

r = Kf + L0 (3) , 6 (3) =L? + K6 {3) , 

the other two basis covectors remaining unchanged. Here 

K=[l- 4l 2 e 2a (r 2 + l 2 )~ l tan 2 ^)]" 1 / 1 , 
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L = 2Kle a (r 2 + l 2 )' 1 ' 2 tan(tf/2), 

so that K 2 — L 2 = 1, the standard Minkowski relation which holds in 
general relativity for an orthonormal basis. Now, 

= -{2K 2 + L 2 )U, X (2)(2) = {K 2 + 2L 2 )U, X {m) = U, 

Y (m) = (2K 2 + L 2 )V, Y {2){2) = -(K 2 + 2L 2 )V, Y m3) = -V, 

the tensors ramaining diagonal and traceless, while neither of the fields 
may be transformed away by these means. This situation shows of 
course more parallels with the electromagnetic field case than the de- 
scription of gravitation using only G, uj, and D (see section 5.2), and 
the general transformation laws for X and G can be obtained in strict 
analogy with electromagnetism (see, e.g., Landau and Lifshitz (1971), 
although in gravitational case the transformation coefficients enter qua- 
dratically, as it is above). 

5.7 Example of the spinning pencil-of- 
light field 

As a pencil of light we understand an infinitesimally thin rectilinear 
energy distribution moving along itself with the velocity of light. Thus 
no electromagnetic field is to be considered, and all expressions will be 
taken outside the singular line of the source. In general the intensity of 
the source may vary along the line, these variations propagating with 
the luminal velocity in one and the same direction. The first (approxi- 
mate) study of such a situation was undertaken by Tolman, Ehrenfest 
and Podolsky (1931) (see also (Tolman 1934)); exact solutions having 
been found by Peres (1959) and Bonnor (1969) (the latter has con- 
sidered also the case of a spinning source, and he first gave a proper 
physical interpretation to the results). Later this subject was revis- 
ited by Mitskievich (1981) and Mitskievich and Kumaradtya (1989) 
who have studied dragging phenomenon in the pencil-of-light field (see 
also (Mitskievich 1990)), as well as additivity properties of such fields. 
From these papers it is clear that the pencil-of-light field possesses both 
quasi-electric and quasi-magnetic properties. Formally it may be con- 
sidered as a gravitational wave field, but, as this is the case in the 
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Einstein-Maxwell pp-wave in section 4.5, this field can be in particular 
a stationary one, thus being definitively no wave at all. 

The simplest description of the spinning pencil-of-light field is con- 
nected with a semi-null tetrad, 



ds 2 = 20<°W 1 ) - 9^9^ - 0( 3 >0( 3 >, (5.7.1) 



where 



N 



9(°)=dv, 9W=du + Fdv + Hd(j>, 

N N 

0(2)= dp, 9^=pd<f). 



(5.7.2) 



Here 



F = k\n(ap), H = H(v), (5.7.3) 



H being an arbitrary function of v (as also k and a may be). The 
solution is a vacuum one (outside the singular line p = 0) so that the 
only non-trivial curvature components 

N N k N H 

i?(0)(2)(0)(2)= - i?(0)(3)(0)(3)= J R(0)(2)(0)(3)= 

are at the same time components of the Weyl tensor. 

For our aims it is better to use an orthonormal tetrad which reads 

0(0) = 2- 1 / 2 [(l + F)dv + du + Hd(l)}, 6®= dp, \ 

9 (i) = 2 -i/2[(i _ F )dv -du- Hd4>\, 9® = pd4>, } {bJA) 

and the non-trivial curvature components are 

fl(0)(2)(0)(2) = #(1)(2)(1)(2) = #(0)(2)(1)(2) = --R(0)(3)(0)(3) = --R(l)(3)(l)(3) 

= --R(0)(3)(l)(3) = k /( 2 p 2 ), 
#(0)(2)(0)(3) = #(1)(2)(1)(3) = fl(0)(2)(l)(3) = %)(2)(0)(3) = ~H/{2p 2 ). 

Thus in the reference frame with r = 

X(2)(2) = -X(3)(3) = F(2)(3) = k/(2p 2 ), \ 

F(2)(2) = -y(3)(3) = -X(2)(3) = ij/(2p 2 ), J ^ ^ 
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or for (anti-)self-dual representation, 



/ \ 

X= 1 ±1 {k±iH)/{2p 2 ) 
V ±i -1 / 



(5.7.6) 



(cf. expression (5.6.9) for the NUT field). Here the numerator is an 
arbitrary complex function of the retarded time v. The field resembles 
that of a charged rectilinear wire with a current in electrodynamics, 
but the gravitational case is much richer, incorporating in particular 
such a property as angular momentum ("polarization") of the source 
(in electrodynamics this would be probably an effect of the Aharonov- 
Bohm type). 

The monad field r = 9^ (5.7.4) corresponds not only to an accel- 
eration, but also to a rotation, 



G = (-k6<?> + H6^)/{2p), oo = -{HO® + k6®)/{4p), (5.7.7) 



which are here orthogonal to each other. We see that in this representa- 
tion, as it was the case for (5.7.5) too, the both linear and angular mo- 
menta of the source (neither of which is transformable away by switch- 
ing to frames in translational and/or rotational motion), contribute to 
the both quasi-electric and quasi-magnetic parts of the total gravita- 
tional field. But in this approach characteristic to section 5.2, the prin- 
ciple of equivalence permits to get rid of rotation or acceleration, and 
moreover, not only locally, but also in a finite region (we speak then cor- 
respondingly on either a normal, or geodesic congruence), so that only 
in the gravitational inhomogeneity field approach, there is guaranteed 
the invariant property of a field to be of quasi-electric, quasi-magnetic, 
or null type, independent of a concrete choice of reference frame. This 
shows once more that for gravitational fields the feature of inhomo- 
geneity plays a crucial role. We may thus choose a non-rotating monad 
congruence, e.g. as r = N(dv + du), N~ 2 = 2(1 — kln(ap)) — (H/p) 2 , 
but this only makes the expressions for X and Y unpleasant, so we 
shall not consider here this subject in more detail. It is however worth 
mentioning that a non-rotating frame can be simultaneously introduced 
in the space-time of a pencil of light only in a limited band of the radial 
coordinate p. This is the so called local stationarity property known 
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also for the ergosphere region of a Kerr black hole (cf. Hawking and 
Ellis (1973)). This limited band can be moved arbitrarily, but it cannot 
be expanded to fill the whole space-time of a pencil of light at once. 

5.8 Gravitational fields of the Godel uni- 
verse 

Finally we come to the Godel space-time (Godel 1949) which has certain 
remarkable properties. We consider its metric in the form 

ds 2 = a 2 [(dt + V2zdx) 2 - z 2 dx 2 - dy 2 - z~ 2 dz\ (5.8.1) 

with the corresponding natural orthonormal basis. It is clear that the 
co-moving tetrad is r = a(dt + \/2zdx) which is rotating but geodesic 
(G = 0, uo = 2~~ 1 l 2 dy) ) while the stationarity yields absence of defor- 
mations (D aj 3 = 0). We assume here that the cosmological constant is 
equal to zero, so that /i — p — (2/ta 2 ) -1 (see, e.g., (Stephani 1990)). 
Then the intrinsic gravitational inhomogeneity field is very simple, 

-2X(i)(i) = X( 2 )(2) = -2X(3)( 3 ) = (6a 2 )" 1 , Y{ a )(p) = 0. (5.8.2) 

Thus there exists no quasi-magnetic gravitational field in the Godel 
universe, although this is rotating (in a sharp contrast with the electro- 
magnetic field in generalizations of this universe where a charged fluid 
does not produce any electric field in its co-moving frame, see section 
4.4). 

As to the gravitational inhomogeneity field sources, they may be 
seen from the quasi-Maxwellian equations (5.4.3) with a substitution 
of (5.4.8)— (5.4.11) and (5.4.16)— (5.4.19): 

X aa > a b° 5 = («/3)av&?, (5-8.3) 

-3u a X aS = k(/i + p)u 5 , (5.8.4) 

^(X$e aSl + Xte aPl ) = 0, (5.8.5) 

while the equation corresponding to (5.4.11), reduces to = 0. Equa- 
tions (5.8.3) — (5.8.5) are satisfied also indeed. We have already ex- 
cluded from them terms being identically equal to zero, this is the 
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reason why (5.8.4) and (5.8.5) are algebraic with respect to X. In fact 
(5.8.4) contains an analogue of div B in Maxwell's equations, and it 
shows how the other terms mutually compensate to permit Y to vanish 
completely. Similarly (5.8.5) corresponds to Maxwell's equation with 
curl B. 



Chapter 6 

Concluding remarks 



In this paper it was of course impossible to sensibly characterize the 
whole variety of use of reference frames in general relativity, unclud- 
ing their new development tendencies. We have deliberately ignored 
here the case of normal congruences (theory of holonomic hypersur- 
faces) which is so important in building the canonical formalism in 
general relativity to the end of gravity quantization. Elementary ideas 
of quantization are discussed in (Mitskievich 1969) and fundamentals 
of the canonical formalism, in (Misner, Thorne and Wheeler 1973, Mit- 
skievich, Yefremov and Nesterov 1985). It is however worth mentioning 
here some new works which could be helpful in understanding the recent 
trends and developments in these areas. 

Much difficulty in realization of canonical formalism in field the- 
ory is due to the problem of equations of constraints (see Mitskievich, 
Yefremov and Nesterov 1985). With help of gauge theory ideas, and 
in particular those of the Yang-Mills theory, Ashtekar (see Ashtekar 
1988a) succeeded in making exceptionally good choice of field canon- 
ical variables and applying them to different types of problems (e.g. 
Ashtekar 1986, Ashtekar 1988b, Ashtekar 1989). One of trends in the 
development of this theory consisted in a study of spaces with self-dual 
Weyl tensor (Koshti and Dadhich 1990) for which a set of general theo- 
rems were obtained (in a more conventional approach, such spaces are 
often called "heavenly spaces" (or 7i-spaces); they were studied and 
applied to solving Einstein's equation, especially by E. Newman, J. 
Plebahski, and D. Finley). The very foundations of the gravitational 
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field theory were also considered, and among them the action principle 
with the Palatini approach (Capovilla, Dell, Jacobson and Mason 1991, 
Floreanini and Percacci 1990) and the covariant action formulation (Ja- 
cobson and Smolin 1988). Sometimes new fruitful developments were 
achieved in the theory of initial data using conformally flat spaces and 
space-time complexification (Wagh and Saraykar 1989). 

The attention was focused in this direction however mainly to the 
canonical formalism in gravitation theory where vacuum constraints 
were brought to a solvable formulation (Robinson and Soteriou 1990) 
and canonical transformations were introduced together with the cor- 
responding generators in the field theory (Dolan 1989). Ashtekar's 
formalism was connected with the tetrad approach (Henneaux, Nelson 
and Schomblond 1989), then a first order tetrad form was introduced 
for description of the action principle and formulation of the canonical 
theory (Kamimura and Fukayama 1990), and elegant results were ob- 
tained in a further use of three-spinors in Ashtekar's formalism (Perjes 
1990). 

Several papers, including the initial works of Ashtekar himself, were 
dedicated to development of methods of 3+1-splitting of space-time 
(Ashtekar 1988b, Ashtekar, Jacobson and Smolin 1988, Wallner 1990). 
On basis of the canonical formalism were also, quite naturally, consid- 
ered the energy problems (e.g., its localizability) (Nester 1991). 

Another branch of these studies was directed to obtaining new so- 
lutions of the gravitational field equations with help of specific choices 
of Ashtekar's canonical variables. Hereby solutions for gravitational 
instantons were found (Samuel 1988), spherically symmetric cosmo- 
logical models (Koshti and Dadhich 1989) as well as a wide class of 
spherically symmetric problems of general relativity studied (Bengts- 
son 1990), while in the latter paper a general discussion of further 
prospects in applications of Ashtekar's formalism was given being of 
considerable independent interest. The formalism was applied to the 
Kantowski-Sachs type solutions with a perfect fluid (Bombelli and Tor- 
rence 1990) as well as Bianchi types cosmological models (Kodama 
1988) (here quantization of these models was also performed using 
Ashtekar's canonical variables). It is significant that Ashtakar's for- 
malism does not in general imply non-degeneracy of the metric tensor, 
so that perfectly new solutions, inadmissible from the standpoint of 
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the traditional theory, enter scope of the study. For example, several 
solutions with massive scalar field were thus obtained without explicit 
use of a metric (Peldan 1990). 

As we see, even a single, though explosion-like evolving in the last 
years trend of application of reference frames, has produced by itself 
a great number of fine results. This is not always the reference frame 
trend proper, of the type which was discussed below in our paper. The 
authors go often beyond the limits of physical reference frames while 
considering either complexified spaces (such as self-dual ones which 
may have both coordinated, tetrad or spinor representations), or null 
congruences, or else. But in all the cases, far-reaching analogies be- 
tween the theories of gravitation and electromagnetism (or Yang-Mills 
theory) emerged clearly, which were often used by the very authors 
for determination of their own research directions and formulation of 
new approaches. This sets our hopes on the fruitfulness of the parallel 
description of electromagnetism and gravitation, however fragmentary 
and concise it could be done here. It is worth repeating that we tried 
to give here a review of a part of the reference frames formalism from 
the side being used, successfully enough, for description of observables 
and interpretation of physical effects which are calculated initially, as 
a rule, in a four-dimensional form without any relation to specific ref- 
erence frames, so that further interpretation efforts are unavoidable. 

Returning to the presentation of reference frames formalism given 
in this paper, we feel it necessary to emphasize that it differs from pre- 
vious approaches mainly in a more consistent use of modern differential 
geometry techniques. At the same time, new and more convenient no- 
tations were introduced which are as near to the conventional ones of 
the old three-dimensional vector calculus as it was possible to achieve. 
It is even possible that the reader could first take our notations for 
those conventional ones and consider them to not be general covari- 
ant. As a matter of fact, they however are covariant under arbitrary 
transformations of coordinates, and the distinction between different 
reference frames is to be seen in the use of a concrete monad r which is 
of course a real four- vector (or covector: the coincident notations of the 
both cannot create any confusion). The advantage of such notations 
(as div or curl) is in greatly facilitated interpretation possibilities when 
one considers equations and effects of relativistic physics in arbitrary 
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reference frames. 
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